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Oscillatory Flow near the Entry of a Circular Tube 


By H. BULENT ATaBExK and CuIeH C. CHANG, Minneapolis, Minnesota, USA?) 


A EAN eats Beh 


i Gi 


1. Introduction 


TA ON 


Steady laminar flow near the entry of a circular tube has been studied by 
“many authors. One of the simplest solutions is an approximate one due to 
Tare (see SLEZKIN [2]?), p. 351). TaRG’s special assumption is that the velocity 
of the convection occurring in the inertia terms remains the same as it was at 
“the entrance section. This assumption removes the non-linearity from the 
equations of motion, and it leads to a solution which is in good agreement with 
experimental results. In this paper we shall make the same assumption for the 
case of unsteady laminar flow. 


5 2. The Equations and Boundary Conditions 


va 


4 It will be appropriate to use cylindrical polar coordinates 7, 0, x with the axis 
of x lying along the axis of the tube, and the origin of the coordinates in the 
entrance section. By symmetry conditions, the flow will be independent of 0. 
The velocity distribution at the entrance section will be uniform provided the 

: connection between the tube and reservoir is shaped like a smooth bell mouth. 

' Thus the entrance velocity will be a periodic function of time only, and it may 

_be taken as 


4 u(x,7,t)|,-9 = Yo f +S" (a coskmt+b,sinkw t)|. (1) 
%, k=l 


By WEIERSTRASS’S second theorem (ACHIESSER [1], p. 32) any continuous 
periodic function of period 2” (in the variable w?#) can be approximated 
“uniformly by a finite trigonometric sum in this form. The solution of the 

"problem corresponding to (1) will therefore enable one to make calculations 
corresponding to any periodic entrance velocity. 

In the absence of external forces, the Navier-Stokes equations are 


Ou » Saree sy (Ss 0?u 1 Ou | 2 
ie ae rele mae 0% “ae Ty OF Ox? }? (2) 
Ou ov Cie. + 1 Op » (Se ee OPu 3 
paige 7 hen eo or Nat © yr a? Get) (3) 


1) Department of Aeronautical Engineering, University of Minnesota. 
2) Numbers in brackets refer to References, page 201. 


od \ oe, tle 3 7 
‘ 
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equation of continuity is 
du th SORA (d 
fe oe ar ( 


where u and v denote the x and 7 components of velocity respectively, and ts 
| 


Equations (2) and (3) are non-linear because of the inertia terms on thei 
left hand sides, and the system of equations is difficult to solve. We can how 
make a number of reasonable approximations that will simplify the analysis. 
In the first place we assume that the transverse velocity v is very small (it is 
actually zero on the walls of the tube and on the axis), so that (3) may be 
replaced by 


0 | 
= 0. (5) 
Moreover, the transverse velocity gradient will be much larger than the longi-— 
tudinal one in the region near the walls where viscosity effects are most 
important, and so 0?u/0x* can be neglected in (2) compared with the othen 
terms on the right hand side. At this stage we have simply made the boundary, 
layer approximations usually adopted in the analysis of steady inlet length flow- 
ty now follow TarG in assuming that the inertia terms u (0u/O0x) + ue 
n (2) can be replaced by 


Uy +3 (a, coskat+ 0b, sink @ a| Se. x 
k=l 


A justification for this substitution is that it is valid without error at x = 0, 
and that were (u, v) differs most from 


E ( + Dy (a, coskw t+ b, sink w 0) ; o 
k=l 


namely where the flow has become fully developed (i. e. no longer influenced 
by the entry), the gradients 0u/0x and Qu/dr are both zero so that the substi- 
tution is again valid without error. Of course between the extreme positions 
there will be some error in the substitution. With the three assumptions we 
have just made the equations to be solved reduce to 


+ thy vy (a, cosk wt + b, sink w ft) = ~~ 5b 42 Zo (r), (6) 


Ou Loo (FO). > : 
‘is Seer Se (7) 
: 
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y 
q 
where p = (x, t). The boundary conditions are: 
- 


uty lee) (a, coskwt+ b, sink wt)|, 
= 8) 


ia @ aoe 4-0) wick WS ieee i. 


tAA), i= 0. lor 7 =a and #« > 0 (9) 


Se Oa foe) eat eee 


4 all values of ¢, where a is the radius of the tube. 


3. The Developed Flow 


_ Developed flow is the Poiseuille flow far downstream where the velocity 
_ profile has become independent of the x-coordinate. In this case, with 0u/dx =0, 
equation (4) can be integrated to give 


. where C is an arbitrary function of time. In view of the boundary condition 
_v=O0Oatr =a, it is seen that C(¢) must be identically zero. Therefore v = 0 in 
_ developed flow. Thus the equation of continuity is satisfied automatically and 
the equations (2), (3) and (4) reduce to 


Ou 1 Op 0?u 1 Ou 
Op om abhavet * a eer: (10) 
Opa 
=o, (11) 
“with the boundary condition: 
w= at r=a. (12) 


In our problem the velocity at x = 0 is prescribed. Since the flow is incom- 
pressible, the rate of flow through any cross section of tube at any time is equal 
_to the rate of flow at x = 0 at the same time. In other words, 


a 


fad a2 
[urdr= 


0 


= c + D" (a, coskwt+ b, sink w ‘ : (13) 
k=1 


By multiplying equation (10) by 7 and then integrating it from zero to a, 
it follows that 


1 op by 2y (ou 
4 mwa D' (a,ksinkwt— b,kcoska@t) + x (F1,..- (14) 
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In writing this expression (13) has been used except that the value of 7 
constant term in it, namely a? 4/2, has been lost by the differentiation with 
respect to 7. 

By eliminating 0)/0x from (10) and (14) we obtain: / | 


n , 2/0 0 a) 
=O ody (te sink wt — bk coskort) — 7 (oe), ate < ; (+5). (15) 


, 
| 
| 


If wu is considered as the real part of the function ~ this equation can be rewrit- 
ten as ; 


= n ' , UL i 279} 0 . 
= =O Hod, (ta, +°b,) b ether — Ae (5), eee y G eae (14) 


and the boundary condition (12) becomes 
w= Osat 7 =a. (17)} 


Equation (16) suggests a solution of the type: 


Hr, = gar) + Saute) oP" (18) 


Substitution of this into (16) gives 


n 


t — OM ; 20 gy. ae y a dq. that 
Pa [seiko — ony (ia, +b) e+ 2% (Ge) — 2 [ee 
2v (dg 0 dp dg 

we a y ae ie) =9- 


Because of the linear independence of the functions e**®!, this equation and the: 
boundary condition (17) yield the following differential equations for the q’s: 


d dy coe dq 
dr (” dy ) ma. ( dy ae ry (19) 
with the boundary condition 


qo(4) = 0, (20). 
and 


adh, 1 dq, ko 2 (dq, kw 3 
dr ‘i y ay nee yp tk ( ‘; mane (? ay +O), (21) 


with the boundary conditions 


Gi@i= OF k= 1,2) a (22) 
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a he solution of (19) with the boundary condition (20) is 


Baraat ae je giatl OW as) 


To determine C, equation (13) can be used, and it gives 


Ls 


bio 


raha hl 


RN 


so that C = 2 uy. Therefore 
y2 
Gol!) = 2 Us (1 — a ; (23) 


which is the well known result for steady Poiseuille flow. 
We shall now consider (21), for which the general solution is 


PO ATER EN RS ne eee EST ee EE ER UIT OLN 


40 = Ay Jo (3 *Vka 7) +B, Y, (** *Vka | 


; 2 dq;. 
+ Ug (4, — 1 Bg) + ae (Fe) | 


koa 


ea) aah 


— 
F 
2 


_where A, and B, are constants to be determined and 


a=a /2 i (25) 


' Since for ry = 0 the function Y, becomes infinite, and the velocity on the axis 
of the tube must be finite, the constant B, is necessarily equal to zero. From 


_ (24) it follows that 

Ms ad . re OL : 

(i), = aaah 2 1 yea) 
and using this in (24) we obtain 


q 
q(t) = A; J (98 Vea e) a ae a I CL ye x) 


3? Vk 


oa 


+ Uy (a, —10,). (26) 


The boundary condition (22) gives 
| 7 pl ti | (aj, — 0 dy) 
jy (8? Vk a 


_and substitution of this into (26) yields: 


s _ py {fo 2 Ve 11a) = Jo (8? Va a) ih 
qi (7) = to (a, — # ) | 7, (i8® ko) | : (27) 
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Therefore the final result is: | 
p 


fs (22k a v/a) — Jy (132 Vk x) | eel | 
Ja (i? Vk x) 


(a, — 1 O,) 


4. The Entry Flow 


First we shall eliminate # as we did in the case of developed flow. We 
multiply (6) by 7 and then integrate with respect to 7 from 0 to a; this gives. 


a 


af. ae 
wf rar + He [t+ Slancosk eo t+ Dy sink od) je [ure 
0 


0 


From (13) we see that 


and so the expression for 0f/0x turns out to be the same as for developed flow, 
namely 
1 Op 


n : 2v (0u 
a = Ue a, (a, k sink wt— bk coskat) + Fr ( =) : 


Or 


Substitution of this expression into (6) gives 


Ou . - Ou 
ae + Mo 1 bes (a, coskmt+ b, sink w ‘ x7 
n 
: 2v (Ou v O Ou 
ASE, ks wt—b. ke pe OE ie i) i eee fuc8 
Uo or z (a, k sink w b, Rk cosk w bt) - ( ), ee hee (r “4 : 


This is a linear non-homogeneous equation and is to be solved subject to the 
boundary conditions (8) and (9). 

Since developed flow is a particular solution of (29), this equation can be 
made homogeneous by writing 


u(x,7, t) = u(x, 7, #) + u(co, 7, t), (30) 
where w(co,7,?) denotes the developed flow, which is given by (28). The. 
substitution of (30) into (29) and the boundary conditions (8) and (9) yield the 
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following differential equations and boundary conditions for the function 1: 


One “ae f po (agcoskont + dysink od] | | 

i (31) 
ss lg a ls 

4 a ee ee U alors | 

and 

: u(x, a, t) =0, (32) 

1,(0, r,t) = Uy i Ds (a, cosk wt + b, sink ost] ~ 2m fs 7 | 

{ | 7 (33) 

Brie oe PallieyR arta) = ¥, lamynah avui 

—u Re | (a, — 60) [AS : = ' hea 

pais Jy (#8 Vi a) | 

4 Let us seek a solution of (31) in the form 

. , ; 

q Ux(x, 7, t) = p(x, t) d(7) , (34) 


with the variable r separated from x and ?¢. Substitution of (34) into (31) 
followed by separation of the functions y and ¢ gives 


j 

a sr + Uo } Dy (a, cosk wt + b, sinrk w #)| ss 4 | 

- ha (35) 
| 3 [2 (H+ EE 6H -—-. | 
The separated equations are 

i Sabtene-2 (9), 0 


» 


"eee 


Y + Ug c iS a cosk wt + b, sink w ‘ ess +yMyp=0. (37) 


k=1 


If we first consider the homogeneous equation (36), we notice that since the 
‘boundary condition (32) implies the homogeneous boundary condition 


4) =9, (38) 


‘the determination of a non-zero solution for (7) is an eigenvalue problem. 
The general solution of (36) which is bounded at 7 = 0 ts 


| b0) = 4 Jl) + 5x (GP), 
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From this equation it follows that 
dp\ _ 
(i) aa nae. 
With this result the general solution becomes 
(39) 


br) = 4 [flar) — 5 Ava). 


The boundary condition (38) gives 
? 
ba) =A [fiAa)— %, faa] =o. 


Therefore, for a non-trivial solution, it is necessary to have 


Jola.a) — 7, fia) =0. 


In view of the recursion formula 


~ Le) = hel) + Jol) 


this last expression is equivalent to 
jo(Aa) =0. (40) 
O (JAHNKE and 


If y; (j = 1, 2,...) denote the roots of the equation J,(y) 
(41) 


Emve [3], p. 168), then the eigenvalues are 
ee) 
4; = ole 


The corresponding complete normalized eigenfunction system is 
V2 2y2 Joly; v/a) : 
= and = 1 — 2) ar; = a ee 42 
do db; a Y; Joly) J ( ) 


This system is orthogonal with respect to the weight function 7. 
. If we denote by y, the solution which corresponds 


Now let us consider (37) 
to the function ¢,, after making use of (25) equation (37) becomes: 
Ow; : : Oy; v 
a a +2, (a, coskmt+ db, sinkw | at - - oy,=0. (43) 
ariables to &, 7 instead of ¢, Be 


If we (temporarily) change the independent v 


Uy 
gato 


where 
é = t ; 1] == % 
mM 
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7 1 
Sots) — (a,sinkwt—b, coskqwt), (44) 
“- k=1 
(43) reduces to 
i Oy, Ve o 


RP Wat bite 


The general solution of this equation is 


vy = fila) e798, 


where /; is an arbitrary function. In other words 


TERN OE 


Bee ily eee (45) 


Data al 


_ Now we have to determine f; (x — up s/w). For x = 0 from (45) we get 


a 


@ 


f,(— 25) = ylt, 0) exo? (46) 


a a el ci | 


e) 
5 


mAY" 


Ug S 1 y2F(s)/x2 
i(- : =v. F(s), o| ergs! 


(a2) 


ay 


where wt = F(s) is the solution of (44) for the variable w ¢. If we replace s in 
(46) by s — x w/t», we obtain 


ife— 8) = nf He 2), 0] oe 229 


(a) Ug 


\ a ad 


een VE 


_and then (45) may be written as 


if x 2 \ ais as Xo ne PE w 
wilt, 2) = we F(s— <2), of ert.) / (47) 
where x «/u#, is replaced by x a?/a R, and R = uy a/v is the Reynolds number. 
Therefore, to complete the solution it is only necessary to determine y;,(t, 0). 
Before we do this, we may note the following facts about the function 

s = s(t) and its inverse mt = F(s). A differentiable function has a single- 
valued inverse in an interval, if its derivative is greater than zero in that 
‘interval. Therefore the inverse function of s = s(m?#) exists and is single- 
valued if 

ds “ F 

a sin k bs 

dot) pane eay 12, sinh i emt 

In view of (1) this condition is equivalent to saying that the rate of flow 
through the tube should be positive for all time. 


-ZAMP X1I/13 
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Integrating this last equation we see that 


wt 


$ =| en. d(w t) + Const. (44a) 
0 


where “(w t) denotes the average velocity. That is, for a given entrance velocity 
(which is equal to the average velocity) s is related to the total amount of flow. 
Therefore once (1) is given, the graph of the function s = s(w #) can be readily 
drawn. The same graph also gives the inverse function m ¢ = F(s). In particular 


6; 


Pal es S-S(wt) 


2 A-xe/ark 


sil ars 


0 40 80 120 160 200 240 280 320 360 
wt ——> 


Figure 1 
Procedure to find F(s — A) for a given w¢. (1) For a given wt, read s; (2) subtract Ad = xo2/a R 
from this value of s; (3) read the corresponding (F(s — A) on w ¢ axis. 


since in all calculations the argument of the inverse function F is equal to 
(s — x «/a R), the value of the constant which appears in (44a) is immaterial 
and may be taken equal to zero. Figure (1) represents the function s = s(o t) 
for the entrance velocity w/z) = 1+ 0-5 cosm¢ and in the legend the process 
of finding the value of F(s — x «2/a R) is explained. 
Now let us consider the determination of y(t, 0). For this we observe that 
the series A 
D vil, t) ;(r) , (48) 
7=0 
is a solution of (31) provided it converges and is term-wise differentiable in the 
variables x,7,¢ as many times as required by (31). This solution must also 
satisfy the boundary condition at x = 0, namely 


u,(0, 7, t) = S"y,(0, t) dj(r) , 
j7=0 
where 1,(0, 7, ) is given by (33). This enables us to determine y,(0, t) Indeed, 
by multiplying both sides of this equation by 7¢,(r) and integrating with 


re ve 
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Tespect to r from zero to a, we find that 

aa 

: ,(0, 2) = | (0, r,t) blr) r dr; 7=0,1,2,.... (49) 
. : 

4 


For the zero subscript 


(0, t) = te Just ¥,t) rdr = 0 


Ee An eer woe 


“since the rate of flow with velocity u, is 0. By substituting (33) into (49), the 
~;(0, t) for subscripts other than zero are found to be: : 


: 2 n shoot 
yee l24. ly PR (a, — 1b) é*°*]| 
4 ) y; Ug a aXe 2 Gea e oly? | 
; (50) 
: = _V2a 7 +O My [a,-cos (Rk w t — O,;) + Oy sin (k w t — Oh 
4 Vj 
where 
; ' Re 
Mi= eae Ook) - 


palit | 


Then the substitution of (42) and (47) into (48) in view of (50), gives 
. 
pi (* r,t) = —4 Hod a ; 


, 


x ( +My, ee F(s— <5) — 0,; ) +b sin (z BIS xa) 0.) (51) 
\ 


vi eat) Jo (vi ¥1 a) 
x exp [Fs — a 0) a ES aa ; 
Now we are able to write down the expression for the velocity of the 


original flow. Indeed, by inserting into (30) the values of u(x, 7, ¢) and w(co,7,t) 
from the equations (28) and (51) we find that 


(%=75,) jae yi “sie = Jali is 2] cites 
Je v ) bm (52) 


gz 4S : yu [a, cos(k(w t — A) — O,;) + % sin (k(w ¢ — A) — 6.) 
j j k= 


ee ee 4 


~ | 
oe t 
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In this equation we have replaced F(s — x «?/a R) by (wt — A). The value. 
of the phase displacement angle A can be determined graphically as explained 
in Figure 1. When wu has been determined, v can be obtained by integrating (7): 


- Ou 0 : ; 
run —f Se rdr=— 5 [urdr. (53) 
0 0 


Since the rate of flow through any cross section, given by (13), is only a function | 
of time, it is clear that the condition v = 0 at ry = a is satisfied if v is calculated 
‘from (53). As this component of the velocity is not going to be used in the 
following, its calculation will not be given here. 
The pressure gradient: The x-component of the pressure gradient is given 
by (14). Let us calculate (0u/0r), and substitute it in that equation. Differen-_ 


tiating (52) with respect to 7, and noting that 2 J,(y;) = y; Jo(y;), we find that: | 


Ou a Uy Li ; J, (8? Veer Airy Ree 
(5), =— = Ret 3) yk (a, — 1 b,) z = a exp E (4 @t+ = I} 


Te he oo n 
= ~ 1+ ASS poe [a, cos (k(w t — A) — 8,,;) 


Eee 


+ b, sin(R(w t — A) — 7) exp(— 7 A) é 


Substitution of this equation in (14) and conversion of the resulting equation 
into dimensionless form give 


1 OD, = at : 


oe O(xja) ~ Re (a,ksinkot—b,kcosko t) 


Blea linge a eee | 
2 Re pa Vk (a, — i dy) 7 @ryEay OF (kot + “all 


| 1 lee) wt 
—-5 i # = Pa 1 hss [a, cos(k (wt — A) — 0,,) 


4- b, sin(k (wt — A) — 7) exp(— ra All. 
32 


Let us note here that, although it is not apparent from equations (52) and 
(54), the velocity and pressure gradient are periodic functions of time with a 


period of 2 a/m. This can easily be proved by showing that A is a periodic 
function of time. 


2 Ba 8 rae 
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5. The Inlet Length 


Eee 


As seen from (52) the velocity u reaches its developed value asymptotically ; 
that is, it requires an infinitely long tube for the flow to become fully developed. 
For steady flow, it is a common practice to define the inlet length as the distance 
between the entrance section and the location where the ratio 


- u(oo, v) — u(x, 7) | 
+3 
; u(co, 7) =o 


“is equal to 0-01. Here u(co, v) denotes developed flow for the steady case. 
Therefore, the entrance length L is given implicitly by the equation 


u(oo, 0) — u(L, 0) 


(08, 0) = OO: 


_ For oscillatory flow the velocities entering in this equation are functions 
of time, and the equation does not lend itself easily to numerical calculation. 
In order to simplify the calculations the following equation will be adopted 
: instead to define the oscillatory inlet length: 


i 


; u(co, v, t) — u(L, 7, t) | 0G (55) 
: t+2a/o | si 
4 all : | 
: or i u(oo, 7, t) dt | 
- t |r=0 
' Since for the present case 
t+2n/w 
(a) 
a / u(co, r,t) dt =2u, 


v 


equation (55) becomes 
u(oo, 0, t) — u(L£, 9,4) _ 9.94 


2 Ug 


Substitution of the values of u(oo, 0, ¢) and u(L, 0, t) from equations (28) and 
: (52) respectively into this equation yields the equation 


7 gee 1+ 37 M,; [a, cos(k(w t — A) — 6,;) + 5, sin (k(t — A) — 4.) 
Es kal 
fi- st Jexn(-% 4) =00 | 
x [1 rem sats oa A) =0-01 


for the calculation A. For a given flow the only unknown in this equation is 
the phase angle A. If we can calculate A from this equation, we can determine 

the inlet length L graphically from Figure 1. Indeed for a certain w¢ and 
corresponding A, the difference A of the ordinates of the curve s = s(o #) 


(56) 
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corresponding to the abscissas w ¢ and wt — A is equal to L «?/a R. Therefore 
we have 


Lett a. (57) 
a 


Now let us consider the problem of solving (56) for A. The left hand side 
of this equation is an infinite series. This series converges uniformly for ¢ > 0, 
and A > 0. Because of the exponential factor, convergence is very fast, and 
therefore we may approximate equation (56) by taking only the first term of 
the series. Doing so we have 


Fhe a) 0 hyen(- 4) 001, ay 


where 


By (@t) =1 + DM iy [a, cos(kwt— 6,,) +d, sin(kw t— 9x1)]- (60) 
k=1 


Thus (59) gives us a relation between £,(w ¢— A) and A. This relation is pre- 
sented in Figure (2) for different values of the parameter «. On the other hand, 


7 


/ 


6 5 


-- 
- 
- 
ae 
-_- 
RE 


1. 


$60 320 ©2280 


240 200 ~—*160 


Figure 2 


A versus w ¢ for different values of the parameter «. The dotted lines in the figure represent the 
extensions of the corresponding solid lines to the left. The values of A for these lines are obtained by 
adding 360° to the corresponding abscissas. 


knowing the entrance velocity and the value of the parameter « we can present 
graphically the function f, versus w ¢ from (60). Now let us plot this curve on 
transparent paper taking the scales of w ¢ and B, equal to the scales of A and 


il i atts ais 
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By of Figure (2). Superpose this curve on Figure (2) so that the point (w #, 0) of 
the 8, versus w ¢-graph coincides with the origin of Figure (2) and the w t-axis 
is directed along the A-axis. Then the coordinate A of the point of intersection 
of the proper A versus §,-curve (the one which corresponds to a given «) with 
_the #, versus w ¢-curve gives us the phase displacement angle for the w¢ under 


consideration. 
4 


ta 


6. A Simple Example 


Let the entrance velocity be 
: U|,9 = Ms (1+ 0-5 coswt) and a=al/? = 4. 


This information is sufficient to determine the flow in the inlet length as a 
function of #%) and the Reynolds number R. 

For this flow a, = 0-5, k = 1, « = 4. By substituting these values into (52), 

velocity profiles can be calculated for different values of x/a R. These calcul- 


| Ba DpBe 


wt=60° 


Figure 3 
Velocity distribution in the inlet length for u| 9 = Mo (1 + 0-5 cos w t) anda=4. 
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Figure 4 
Velocity distribution in the developed flow for | x-0 = % (1 +05 cosw?t) anda=4. 
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[Steady inlet length 
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wt —> 
Figur 5 


Inlet length versus w t for 7) Pie = Uy (1+ 0-5 coswt) anda=4. 


ations have been carried out and the results are presented in Figure (3). 
Developed flow comes out of these calculations as a side product and it is 
shown in Figure (4) for different values of w t. Figure (5) gives a presentation 
of L/a R versus w t for the same flow. 
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4 Sommaire 
“4 
On étudie l’écoulement laminaire oscillatoire d’un fluide incompressible visqueux 
‘prés de lentrée d’un tube circulaire sur la base des €quations linéairisées Navier- 
Stokes. On suppose que la vitesse du fluide soit uniforme a l’entrée et qu’elle soit 
‘une fonction de temps connue périodique arbitraire. On obtient des formules pour 
les composantes longitudinales de la vitesse et le gradient de pression et on calcule 
‘la longueur d’entrée. 
q Zusammenfassung 
Die periodische laminare Str6mung einer zahen inkompressiblen Fliissigkeit 
‘in der Eintrittspartie eines zylindrischen Rohres wird auf Grund der linearisierten 
Gleichungen von NavierR und Sroxes berechnet. Es wird dabei angenommen, 
‘dass im Eintrittsquerschnitt in jedem Zeitpunkt eine iiber den Querschnitt kon- 
stante Geschwindigkeit herrsche, die jedoch eine Funktion der Zeit ist. Formeln 
fiir den Druckgradienten und die Geschwindigkeitskomponenten in Richtung langs 
‘der Rohrachse werden erhalten. Ebenso wird die Einlaufstrecke berechnet. 


(Received: September 19, 1960.) 
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Ein mechanisches Rechengerat zur Behandlung des ebenen 
Dreikérperproblems 
Von P. Curtt, Ziirich 


Kraft und Weg 


Wir betrachten zunachst nur einen Massenpunkt P,, unter dem Einfluss 
einer Kraft, die standig nach dem festen Zentrum M zeigt (Zentralbewegung, 


Figur 1) und deren Grésse eine Funktion der Entfernung a ist (a4 = MP. xy) 
Es kann dies zum Beispiel die Kraft eines gespannten Gummifadens sein 
(elastische Krafte) oder die Schwerkraft, die stets durch das Erdzentrum geht. 

Handelt es sich speziell um eine Gravitationskraft, so tibt die Masse M auf 
P., eme Anziehung aus, deren Beschleunigung nach NEwron 


b= — pe oF (k = Gravitationskonstante) 


ist. Der Index von by bezieht sich auf diejenige Masse M, welche die Be- 
schleunigung austibt; angesetzt wird der Pfeil am Punkt P,,, welcher die 
Beschleunigung erleidet. 


Figur 1 


Integrieranlage K, Aufzeichnung der Bahn 


Mit den ee von Dy, eeenen sich die Bewegungsgleichungen fiir 
den Punkt P,,: 
dx ss ee BN Li) eee a ME a 


dt® e a ea? (1) 
d*y M M y 
Gi ge is sin a= —Rk7> =. (2) 
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Durch zweimalige Integration nach der Zeit ergeben sich daraus die Wege s: 


[oa=, fode=s. 


‘Wir betrachten zuniachst b, und b, als gegeben und integrieren mechanisch wie 
folgt: 

a) Wir setzen b, um in die Drehung einer Welle W (Analogierechner); diese 
fihren wir zum Integrator j,, der — durch den Zeitmotor unter dem Tisch 
angetrieben — die Integration di b,: dt =v, besorgt. Demnach tritt rechts v, 
aus dem Integrator J, heraus und eeee in if hinein, wo die zweite Integration 
i, v,° dt erfolgt, so dass rechts dx aus J, heraustritt. Analog liefert der Ein- 
gangswert 0, die Weg-komponente dy am Ende rechts des Kastens K,. 

mesa ausgedriickt: Wenn wir links die Beschleunigungskomponen- 
ten b, und 6, in den Kasten K, hineingeben, dann kommen rechts die Weg- 
Be caponeriten dx und dy heraus. 

Das Detail der Integration iibergehen wir, die Kreise kénnen sowohl 
Scheibenintegratoren als auch Kugelintegratoren darstellen, wie solche in 
letzter Zeit von der Firma Amsler bedeutend verbessert worden sind. Verstar- 
ker sind weggelassen, der Zeitmotor, der die vier Integratoren antreibt, befindet 
sich, wie schon erwadhnt, unter dem Tisch. 

Ue sind Ubersetzungsgetriebe fiir den Wechsel der MaBstabe und Verzer- 
tungen, K Kupplungen. 

b) Die x-Welle endigt in einer Spindel Sf, die den Laufer L nach rechts 
oder links bewegt. 

Am Ende der y-Welle ist ein gezahnter Zylinder Z, welcher den Stab St 
(dessen untere Flache gezahnt ist) entsprechend y in den Fihrungen F des 
Laufers L verschiebt. Die Spitze des Stabes St steht momentan auf P,,, si 
beschreibt die gestrichelte Bahnkurve. 


Prinzip der mechanischen Lésung 


Die vorstehend beschriebene Apparatur zur mechanischen Integration ist 
bekannt und kein Problem. Schwieriger ist es, ihre Eingangswerte b, und 0, zu 
bilden, da diese vom Ort P,, abhangen, also von Endwerten x und y, die zu 
ermitteln wir erst im Begriffe sind. 

Es kommt hinzu, dass wir spadter einen zweiten Kérper mit ins Spiel 
bringen werden, so dass dann vier Differentialgleichungen simultan tntegriert 
werden miissen. 

Wir betreten deshalb den Weg einer mechanischen Losung, da bei einer sol- 
chen das Prinzip der Riicksteuerung (im Gegensatz zu den graphischen und 
rechnerischen Methoden) uns erlaubt, Endwerte spontan wieder als Ausgangs- 
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werte in einen Kreislauf laufend einzusetzen, und dies mit betrachtlichen 
Zeitgewinn. 

1. Wir fiihren in Figur 2 von P,, aus iiber M einen Faden auf den (nach 
geteilten) Umfang eines Funktionszylinders A. Dessen Mantellinien (Ordinaten 
geben uns dann laufend die Grésse des variablen Wertes b = F (a) an. Wi 
tasten sie bei A ab, um sie durch eine Kette Ke von unten in die Rille R de 
Richtungsscheibe S tiberzufiihren. Auf der Unterseite der Scheibe S lieger 


Figur 2 


Prinzip der mechanischen Lésung. 


parallel zu R zwei Zylinderfedern, die sich in den Punkten I und II abstiitzen 
und den Endpunkt der Kette nach aussen driicken. 

2. Die Richtung von a (Linie P.M ) nehmen wir mit dem Tastarm T ab. 
Ein Motor 1 dreht die Wellen w und ihre Schnecken Sn und damit die Rich- 
tungsscheibe S sowie den Tastarm 7. Beriihrt der letztere den Metallfaden 
P.M , so schaltet ein Relais den Strom zum Motor 1 ab; so wird auf det 


Richtungsscheibe S die Beschleunigung 0 stets zu J »yM gleichlaufend gehalten 
(Kreisnachlaufwerk). Wir haben damit } nach Grésse und Richtung. 


3. Schliesslich zerlegen wir den Vektor 6 auf der Richtungsscheibe S in 
seine Komponenten b, und b,: 

Die Spindel S#, angetrieben durch den Motor 2, verschiebt den Tastarm 7% 
nach links. Erfolgt Kontakt, so schaltet ein Relais den Strom des Motors 2 ab. 
Die Zahl der Drehungen von S¢ ist ein Mass fiir b,; sie wird der Integrieranlage 
K, durch eine Kette Ke zugefiihrt. 
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Ahnlich wird die Komponente b, abgetastet; wir wollen aber damit die 
‘kleine Figur nicht iiberlasten. 


4 4. Schalten wir zuletzt den Zeitmotor ein, so zeichnet das obere Ende der 
‘Stange St automatisch die Bahnlinie auf. Der ganze Ablauf der Apparatur 
wird emnzig und allein durch das Newtonsche Gesetz regiert; wir brauchen nicht 
Zu wissen, dass die Bahn eine Ellipse, der Hodograph ein Kreis ist und M der 
Brennpunkt. 

_ Vorerst allerdings ist der Apparat keine Notwendigkeit, da wir den Bewe- 
“‘gungsvorgang rechnerisch durchaus beherrschen. 

_ Die Anforderungen an die Genauigkeit sind betrachtlich: Das Ende der 
Bahn muss genau und tangentiel in den Anfangspunkt einlaufen, um den 
‘Umlauf auf der aufgezeichneten geschlossenen ersten Umfahrung zu wieder- 
olen (spater, beim Dreikérperproblem, sind die Bahnen nicht immer ge- 
schlossen). 


q Zwei Kérper M und uw 

Grosse der Beschleunigung b (Figur 3): 

_ 1. Die Krafte K, mit denen zwei Massen M und yp sich anziehen (Figur 3a) 
sind gleich, namlich 


Zz ees ees 


(Newton) +) 


2. Dagegen sind die Beschleunigungen ungleich Figur 3b ; die grosse Masse 
M iibt auf jedes Teilchen der kleinen eine grosse Beschleunigung aus, die kleine 
Masse ys vermag den grossen Kérper M nur wenig zu beschleunigen. 


u f ee 
ws by-K- Fe rd ee : 
V4 


Figur 3 
a Krafte b Beschleunigungen c Zusatzbeschleunigung 6, 


1) Aus dem Getriebeplan ist ersichtlich, wie unser Analogiegerat ein tibersichtliches Modell des 
Bewegungsvorganges ist: wir sehen auf dem Reissbrett die drei Massen M, uw, m ihre Bahnen ziehen, 
ferner auf dem Brett der sechs Richtungsscheiben S laufend Grésse und Richtung der sechs 
Beschleunigungskomponenten sich andern. 

Nun ist das Dreieck M yu m allerdings eine Darstellung von begrenzter Grosse, welche hin- 
sichtlich grésserer Genauigkeit keine bedeutende Steigerung erlaubt. Um eine solche zu erreichen, 
miissen wir die Dreieckseiten a, b, c und die Richtungen «, B, y ebenfalls als Drehungen von Wellen 
(Analogierechner) darstellen ane auf engem Raum sozusagen aufspulen. Das fiihrt auf eine sehr 
umfangreiche Apparatur. 
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Ist M die beschleunigende Masse, so iibt diese eine Beschleunigung aus: 


bu = bie ebenso ... 0, = hp 


3. Beide Massen M und y sind in Bewegung. Wollen wir M zum Stillstand) 
bringen (Figur 3c), so miissen wir in M eine gleichgrosse, in der Richtung: 
entgegengesetzte Zusatzbeschleunigung b, anbringen, und (damit das Systen 
der beiden Massen in sich ungestort bleibt) ebenso in w (Figur 3c). . 
Betrachten wir also die Bahn von yu in bezug (relativ) auf die feststehende 
Masse M (im Koordinatennullpunkt N), so haben wir ihr nach Figur 3c eine 
zentrale Relativbeschleunigung . 


Bret = Ox oF b, =k ae fe F(a) 


zuzuteilen; sind die Massen gegeben, so legen wir 
Dra = F(a) 


um den Funktionszylinder A von Figur 2. Die vom Apparat Figur 2 aufgezeich 
nete Kurve ist dann die Relativbahn von yw, betrachtet vom festen Stando 
M aus. 


Drei Koérper 


Die Beschleunigungen: 

Die drei Massen sind M, uw, m; iiber ihr Gréssenverhaltnis wird keine ver 
einfachende Annahme gemacht. 

Die Seitenlangen des Dreiecks sind a, 8, c. 

Wir halten jetzt die Masse M fest und betrachten von ihr aus die Bewegung: 
der Masse w und nachher die Bewegung der dritten Masse m. 


Figur 4 


Bewegung der Masse wu. 
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Fiir die Masse wu gilt nach Figur 3c: 


Peas =i 


a2 


i Nun tritt aber die dritte Masse m als Stérmasse hinzu und erteilt der Masse le 
- eine Beschleunigung 


bo = (kG) 


- und der Masse M eine Beschleunigung 
n 
bo =(e 7) I 


_ Da wir aber M festhalten wollen, muss dort zu dy; in entgegengesetzter Rich- 
_ tung eine gleichgrosse Zusatzbeschleunigung III (punktiert) angesetzt werden 
_ und ebenso (damit das System der drei Massen in sich nicht gestért wird) 
_ bei w und bei m. 

Fiir jede der drei Beschleunigungen bendtigen wir einen Funktionszylinder 
(I, II, III); andern die Massen, so sind die Funktionslinien neu aufzuzeichnen. 


Storglieder 


2. 

a 
‘Tar (x4) I (K-) iia 
mt a mm ” a 


Storglieder 


Figur 5 
Ubertragung der Distanzen a, b, c (Fadenlangen) auf den Umfang der Funktionszylinder. 


Die drei Zylinder werden (siehe Figur 5 und auch Getriebeplan) durch die 
Metallfaden a, b,c gedreht. Ein Faden, der in m befestigt ist, ist iiber wu und 
M auf das Differential D gefiihrt, wo er das eine Sternrad gemass (a + 0) 
dreht; das andere Sternrad wird vom Faden a gedreht, so dass uns die Achse 
des Planetenrades im Differential die halbe Differenz [(a + b) — a]/2 = 6/2 gibt 
und hierauf die Ubersetzung 2.1 die verlangte Drehung 6 des Zylinders II. 
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Fiir die Masse m und w als Stérungsmasse sind nach Figur 6 die Beschleunigun- | 
gen im Punkt m: 


bra = (2 an) IV, by=(k5) V, b,=(k4) VI. 


: 


Figur 6 


Bewegung der Masse m. 


Die Funktionszylinder sind IV, V, VI in Figur 5. Dabei sehen wir, dass I und 
VI vom Faden a gedreht werden, II und V vom Faden b, III und IV vom 
Faden c, weshalb je zwei auf dieselbe Achse gesetzt sind. 

Von den Zylindern sind die Beschleunigungen durch eine Kette auf die 
Richtungsscheiben S zu libertragen (vergleiche die Prinzipskizze Figur 2). 


Die Ubertragung der Winkel a, B,y auf die Richtungsscheibe S zeigt 
Figur 7. Jeder Motor M treibt zwei Scheiben an und einen Tastarm T. Beriihrt | 


Figur 7 
Ubertragung der Richtungen «, B, y auf die Richtungsscheiben S. Diese z 
Richtung der Beschleunigungen I bis VI. 


eigen laufend Grésse und 


PNeA 
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+ a < 


die Spitze von T den Metallfaden, so schaltet ein Relais den Strom zum 


a 


“Motor ab. 

4 Die Rillen der Scheiben S sind immer den Dreiecksseiten a, b, c parallel. 
_ Mit den Hilfsfiguren 2, 5,7 (zu Figur 9) diirfte nun der Getriebeplan am 
_ Schluss verstandlich sein. 

3 Die Bewegungsgleichungen?) 

_kénnen jetzt aus den Figuren 4 und 6 abgelesen werden. 

_ Fiir die Masse w, deren Relativbewegung wir von M aus verfolgen wollen, 
‘s ind die Komponenten in der X- und Y-Richtung mit m als Stérmasse nach 
Figur 4, 

4 dx M+ 

a 9 i 

: Ee Apes —(k = ) cosa + (k Fr) 008 B — (2 ar) COSY, (1) 
4 

“3 

A a : : : 

; | oe |, = 5 = — (1) sina — (1) sin 8 + (II) sin y. (2) 


Die ersten Summanden dieser beiden Gleichungen bezeichnet man als die 
_Kepler-Gleichung, den zweiten und dritten zusammen als die St6érfunktion. 
‘Die Klammerwerte finden wir auf den Funktionszylindern I, II, III, Figur 5. 
3 Fiir die Masse m mit Stérmasse yu sind die beiden entsprechenden Glei- 
-chungen aus Figur 6 (f nach Figur 4): 


ax M+ 
dt? Cm ire (2 "7 ™) cos y — (ef) eos B — (2 “) cosa, (3) 
2. 
Se |, + by = + (IV) sin y + (V) sin B — (VI) sin a. (4) 


‘Diese vier Gleichungen (1) bis (4), die langst bekannt sind, werden also durch 
‘unseren Apparat kérperlich dargestellt und — wenn wir den Zeitmotor laufen 
lassen — als System simultan integriert. Als Resultat erhalten wir die Relativ- 
bahnen der beiden K6rper mu und m in bezug auf die feststehende Masse M. 


Betspiel (Spezialfall) : 

Die drei Massen M, m, « bedeuten der Reihe nach Erde, Mond, Satellit. 

a) Fir den Mond m gilt Gleichung (3); w als St6rmasse ist jetzt verschwin- 
dend klein, so dass die zwei Stérglieder mit w in der Gleichung (3) wegfallen. 
Die Zylinder V und VI drehen sich zwar (weil die Fadenlangen a, }, c laufend 
sich dandern); da aber die Ordinaten der Funktionslinien auf den Zylindern 
Null sind, bleiben die Abtaster A stehen. So bleibt in Gleichung (3) nur der 
erste Summand. 


2) Vergleiche Prof. Ww. Scuaus, Die himmelsmechanischen Grundlagen der Rawmfahrt, Gl. (78), 
S. 63, im Buche « Raumfahrtforschung» von H. GARTMANN (Oldenbourg-Verlag, Miinchen 1952). 
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b) Beim Sateliten w dagegen gelten in Gleichung (1) alle drei Summanden, 
nur im ersten ist ~ wegzustreichen. Alle drei Zylinder I, II, III arbeiten. :| 
Der Satellit ist senkrecht gestartet in programmgesteuerter Aufstiegbahn. 
in das Perizentrum P der Fahrt-Ellipse eingefahren worden. Dort hat er jetzt 
einen Abstand g, vom Erdzentrum und eine Abgangsgeschwindigkeit v», deren) 
Grésse zwischen der Kreis- und der Parabelgeschwindigkeit liegt und so gross: 
sein muss, dass das Apogdum jenseits der Mondbahn sich befindet. 
Mit v9 Q ist die storungsfreie Ellipse bestimmt; aber die Pfeile 1, 2, 3) 
deuten an, wie der Mond die Tendenz hat, den Satelliten aus dieser Bahn zu 
zerren (Stérungsmasse m des Mondes in der Satelliten-Gleichung (1)). 
Wir kuppeln jetzt aus und stellen am untern Rand rechts des Getriebe-: 
planes ein: 
Am Zeit-Indikator des Zeitmotors: ¢ = 0. 
Mit Kurbel G, an der Spindel Sp: — x = gy und v, = 0 am Integrator J. at 
Mit KurbelG, am gezahnten Zylinder Z: + y = 0 undv, = vam Integrator J,.. 
Ebenso ist am obern Rand des Getriebeplanes, exakt zur selben Zeit ¢ = 0,. 


die Konstellation des Mondes m nach Ort (x, y) und Geschwindigkeit (v,, Vy) 
einzustellen. 


Figur 8 
M Erde m Mond w Satellit. 


d) Legen wir speziell die X-Achse in die Verbindungslinie der drei Massen. 
M,m,wu (Figur 8), dann haben wir in Figur 4 fiir den Satelliten f die drei 


Beschleunigungen I, II, III nach links in die Horizontale zu drehen; ihre 
Summe ist dann 


S=l+ +I] —ette pp ™@ opm 


a® pe OER 


Die Richtungsrillen R aller Richtungsscheiben S liegen jetzt waagrecht und 
alle drei Faden a, b, c liegen in der X-Achse. Bezeichnen wir noch (Figur 8): 


(TK 


? 
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_ so folgt: M+u m m 
=|b|,=—k = =e 
A may | le Re o ge! 


das) 


- wie bei SCHAUB, S. 86, Gleichung (122), Fall 2, wenn je= 0. 

Die sekundaren Einfliisse (Sonne, Erdabplattung) sind auf zusatzliche Be- 

 schleunigungen ¢, ¢, zuriickzufiihren und diese durch Differenziale auf den 
beiden Wellen fiir 6, und 6, zu addieren. 


ee 


COENEN nas aes 


: : 
4 

3 Vv Vv Wl 

4 lox —-{k F) cos 7-(Kjz) cos B-(K-57) cos a i 
’ 

} Vv Vu 

i ee i 

; ayn 

4 

4 

a I~x,¥ a 


rT ! 
Ds 

'y 

| 

=] Sri 


ee I 
26, 
oS by 
Figur 9 
Getriebeplan. 
Summary 


The project concerns the problem and principle of solving the three body 
problem by known mechanical means: function cylinders with a photoelectric 
follower, integrators, and adding machinery. The equations of motion each contain 
three terms (including the perturbation function); these are represented indivi- 
dually on the function cylinders, resolved along the axes x and y, and plotted in 
the working plane. 

One of the three bodies is held firm, and the apparatus then records auto- 
matically the relative paths of the other two. 


(Eingegangen: 23. August 1960.) 
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On the Centre of Flexure : 


By Emeric Deutscu, Bucharest, Rumania?) 


1. In a paper printed in 1950, BERMAN [1]2) has determined the coordinates | 
of the centre of flexure of a prismatic beam when only the boundary-values of | 
the torsion function are known. In this paper we shall arrive at this result in a | 
much simpler way and without taking the axes to be the principal axes of the 
cross-section. We shall apply the result thus obtained to the case when the 
section of the beam can be conformally mapped on the unit circle. 


2. Let a cantilever beam of uniform, simply-connected cross-section have 
one end (z= 0) fixed and the other end (z = J) loaded by a force (P,, B, OF 
lying in the plane z = /. The z-axis is taken along the central line of the beam, 
while the x- and y-axes are any orthogonal axes intersecting at the centroid of 
the end z = 0. 

We denote the boundary of the section z = 0 by J” and the region interior — 
to I’ by R. We introduce the notation 


EAM ii “am yr dx dy . (1) 
J | 


We have, since the origin is the centroid of R, 
Ii => Lo => 0 e 


It is known [2] that the flexure problem is reducible to the problem of 
finding three functions g, y, and q,, harmonic within the region R, whose 
normal derivatives are prescribed on the boundary: 


ee = y cos(n, x) — x cos(n, y) on 22 7 (2) 
Parakey) ghey 21¢ YR 

“ y?] cos(m, x) on : (3) 
Ope ei) 2 2 

Fn lA +») 2 — 9» 28] cos(m, y) . on I. (4) 


Here m denotes the exterior normal to the boundary J’, y is Potsson’s ratio and 
¢y is the torsion function. 


1) Professional School », Klement Gottwald”. 
*) Number in brackets refer to References, page 218. 
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The coordinates (x,, y,) of the centre of flexure are given by 


%.= J (La S2— TyiSy), Y= TS (111 S2— Iy,2 Si) » (5) 
where 

J7*=2(1 +») (Io,2 tog 123), (6) 
and ; 


ol [= Ge —y Pt +9) Py — vy] dx dy, 
(7) 


oe pints. & —y  — (L49) x+y] dedy, 


We try to express the quantities S, and S, with the aid of the torsion 
function m. We denote 


Dy = ff (« Ze — y Pe) de dy (k= 1,2). (8) 


; _ According to GREEN’S formula, we have 


OS SN Ne ONS) eo 


ie “IES te = 2 | dx dy =~ — o% (x dx + y dy) 
L 


= — % Ly cos(n, x) — x cos (n, y)] ds, 
fp 


_ or, making use of (2), 


Dari OP de ea 22). (9) 


Now, if we make use of the identity 


Ove 0d 
fo 5 as= hp ¥ or as, (A) 
Cc Cc 


where ® and Y are two functions harmonic within a closed region, bounded by 
a simple — closed contour C, the above formula (9) becomes 


D.=-$ 9 ds (k=-1,2), (10) 
Hw 
or, with the aid of (3) and (4) 


ee f p [(1 +») «2 —» y)] cos (n, x) ds (11) 
HAS 


and 
Deep wlth + vy) y2— v x2] cos(n, y) ds. (12) 
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But we have 


> (x? — y*) wm cos(n, x) ds 
r 


OC ae Py 
Sa ge Oo) re eee 
= (4 -« 2) SE as 2G yxy cos(n y) ds 
m7 3 y on j 4 
P 


(= — x9) (x dx + y-dy) +2 f px y cos(n, 9) ds 
r 


Il 


314 — Tos + 2 px cos(n, y) ds 
? i 


where we have successively made use of formula (A), boundary condition (2), — 
GREEN’Ss formula and notation (1). 
Thus, (11) becomes 


D, =» (Io,3 — 3 Ig1) — f y [x2 cos (n, x) + 2» xy cos(n, y)] ds; 
F 
similarly we find 


D, = v (3 Ih. — I3,9) — f y [2 xy cos(n, x) + y? cos(n, y)] ds. 
P 


With the aid of these expressions for D, and D,, relations (7) become 


S; = (1-27) In, — f @ [a® cos (n, x) -+2yxycos(n, y)] ds, | 
He 


(13) 
S,=— (1-2) 1,,.—- f y [2x ycos(n, x) + y? cos(n, x)] ds. | 
vi 


Unlike the formulae (7), relations (13) help us to find the quantities S, and 
S. (and therefore the coordinates x,, y,), without our necessarily knowing the 
functions ¢y, and g, but only the values of the function g on the boundary J”. 

If the axes are taken to be the principal axes of the region R, we have 
Z,,, = 0 and thus formulae (5), with the aid of (13), become 


2(1 +) Ipex, = —(1 — 29) Iy2—- f p2v xy cos(n, x) + y® cos(n, y)] ds , 
tr 


2 (1+) logy, = —(1—2) In, 4+ f @ [2* cos (n, a) +2xycos(n, y)] ds. 
ve 


These formulae coincide with those found by BERMAN by a more difficult 
method. 


It can be shown that formulae (13) may be put in the form 
4 S,=(1—2»)I,,—-ImQ, S,=—(1—2»)I,,.—ReQ, (14) 
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q Q= 5 f (Me) + /@) [0 — 2») Baz + (+29) Pde—227az. (15) 
: 


; z=x+tiy and f(z) =9(x,y) +iylz,y) 


Es the complex torsion function. 
_ Formulae (5) and (14) reduce the problem of determining the centre of 
flexure to the calculus of the integral (15). 


3. Let 
a(¢) , (16) 


ibe the function which maps the region R on the unit circle | €| <1 in the 
¢-plane, and let 


; folC) = fleo(0) 


denote the complex torsion function expressed in terms of the variable ¢. 
_ It is known that 


hi) = x5 § MSE ao, (17) 


Y 


where y is the boundary of the circle | ¢ | = 1 and ao = e’”. 
Formula (15) becomes 


O= 5 ¢ [flo +f (2)] | 
x [« ay ») w2(a) deo(a) + (1 + 2) o(-) don(o) — 2 w(a) o(=) do )| | 


Formulae (17) and (18) yield the value of Q if the mapping function z=a(¢) 
is known. 


4. Suppose the mapping function (¢) is written in the form 
2 = w(C) = itt, Ee, (19) 


a, (n = 0, 1, 2, ...) being complex coefficients. Let us find in this case a formal 
expression for Q, in terms of the coefficients 4@,,. 


(18) 
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It is known that 


fC) =i D7 8, 
n=0 


where 
+00 3 
b, in Dia a, e 
r=0 
It is found, after some algebraic manipulation, that 


= 3%, 0", w(a) dw(o =i 56, o" do, 


aT) nn paw ct em | ave 
BH {) dolo) = 1 JB, ord, 20) Wl 5) dol) = —1 D/C, 0" 48, 
where 
a, =), 4,-,% (n=0,1,2,...), 
p=0 
By = P Ay Xp (n= 1,2 Rs 
p=1 
+00 
B,=)D/ (b+ 2) G@,1, (*=0,1,2,...), 
p=0 
+00 
Bove SIDR ee rete 
p=1 
C= Pa a, a0, been, | 
p=1 
Ca, =) (2 1) a, 4,4 (eee | 
p=0 
According to (22), we have 


f [folo) + fo( 5) %(e) deo(a) = 2 370, By 


Y n=l 


f [alo + if )] (5) dole) = 22 3, B_, ~ 5, B) 


2 p [falo) + f(-+)] (0) a( +) do(+) 20 3°, C_, 5, 6,), 


20) 


(23) 


(24) 


(25) 


(26) 


ae 


me 
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See a 
Dn, 


since o = e’® and the integral of every term involving o” (n + 0) vanishes. 
Thus (18) becomes 


e 


Q= 7) (0 29) b, By - (1 +27) (b, Bo b, BS) “iF SOs Om 
n=1 


aa i kL 
om 
la} 


28s) a 
Q= = J", By — b, An + , An — 2 (0, B, + 5, B_, — 8, B,)], (27) 
a=) 


“where 

5 Av=B,+C,=)) (2p+7) a a,., (1 =0,1,2,...), 

4 =0 

vonpe (28) 
‘ eee Cee, (P41) Cp np b= 1, 2h idee 

d p=0 


Obviously, formula (27) may be directly applied whenever the mapping 
function w(¢) is a polynomial or whenever the mapping can be approximated 
with sufficient accuracy by a polynomial. 


5. As an illustration of the foregoing procedure, we consider a beam whose 
_ cross section is a cardioid. When the origin of the axes is at the centroid of the 
section, the parametric equations of the cardioid are 


a 


x=2+6cost—3cos2t, y=6sint—3sin2?¢. 


| The analytic function which maps the interior of the cardioid on the unit 
peircle |¢ | = 1 is 
| Z=oC)=2+66—37¢. 


The non-vanishing coefficients a,, b,, &,, B,, A, A_,, B,, B_, are given in 
the following table: 


~ 
¥ 
| 
r 
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They were found after some straightforward computations with the aid of 
(21), (23), (24), (25) and (28). 
We have now 


x6, Bn = OB, cy b, A,) = 9/2 , » (bn a. ak: b, Dies, _ b, B,,) = —2268 , ; 
n=1 neo 


and therefore from (27) we have 
O=8ia (3-147). 


According to (1) we find 


gs ets a, Io,=]=0,, Ig= =o, Le 


oo 
I 
= 
~ 
o 
nr 


ye 


17012 
a=, N= 0h 


From (6) and (14) we have now 


fF 221701 & V4) dg gs StH 0) Spiel a Gee 


and now from (5) 


ad 2(3+ 47) 2 
a Gers mCi Ce 
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Résumé 


Dans une publication publiée en 1950, BERMAN a déterminé les coordonnées du. 
centre de flexion d’une barre prismatique lorsqu’on connait seulement les valeurs 
de la fonction de torsion sur le contour. 

Dans le travail présenté ici on retrouve ce résultat par une voie beaucoup plus 
simple, et ceci sans choisir les axes principaux comme axes de coordonnées, 


On applique le résultat trouvé au cas ot la section de la barre peut étre repré- 
sentée d’une maniére conforme sur le cercle unité. 


(Received: October 3, 1960.) 
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_Hydromagnetic Flow over a Plane Wall with Uniform Suction 
: By TsunEHIKO Kaxuranl, Kyoto, Japan?) 
; 1. Introduction 
4 This paper deals with the flow of an incompressible, viscous and electrically 
_ conducting fluid over a non-conducting plane wall with uniform suction in the 
_ presence of a transverse magnetic field. 
, Assuming that all physical quantities depend only on the distance y from 
_ the wall, we can obtain an exact solution of the hydromagnetic equations, i. e. 
the simultaneous equations of modified Navier-Stokes and MAxWELL’s 
Z _ equations. The solution thus obtained is valid or fails to be valid according as 
the magnetic pressure number S, based upon the suction velocity and the 
intensity of the applied magnetic field, is smaller than unity or not. The velo- 
_ city distribution reverts, of course, to the ordinary hydrodynamical one when 
_ S is equal to zero and/or the magnetic Prandtl number P,, is zero. 
It should be noted here that no steady solution can be obtained for the case 
of uniform imjection instead of uniform suction being distributed upon the wall. 
From the results of the calculations, it is found that both the shear-stress 
and the heat-transfer at the wall are not affected by the hydromagnetic effect, 
and that the displacement thickness of the boundary layer increases with S 
but does not depend on P,, unless P,, is zero. It is quite surprising and somewhat 
paradoxical that the displacement thickness does not depend on P,,. This point 
will be discussed in the main text. 
Recently, the same problem was treated by Gupta [1}?) who concluded 
that the steady solution can be obtained for both cases of suction and injection 
when S > 1, and no solution can be obtained when S < 1, in contradiction to 
the present result. It seems, however, that his boundary condition for the 
magnetic field at the wall does not reflect the fact that the magnetic field is 
applied transversely. In fact, his solution gives a finite value of H, at the wall, 
H, being the component of the magnetic field along the wall. On the other 
hand, since the wall is assumed to be electrically non-conducting as in the 
present paper, no electric current flows in the wall and thus the magnetic field 
in the wall should be curl-free, which leads to H, being a constant in the wall. 
This constant must be identical with the value of H, at the wall by virtue of 


1) Department of Physics, Faculty of Science, University of Kyoto. 
2) Numbers in brackets refer to References, page 230. 
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the continuity condition of the tangential component of the magnetic field at 
the boundary. This means that the magnetic line of force in the wall is not 
transverse but oblique to the wall. Thus we may say that his solution is an 
isolated one in a sense that it is valid only under a particular magnetic field 
whose x-wise component in the wall is just identical with the boundary value 
given by his solution. Mathematically speaking, since the highest order deri- 
vative of H, with respect to y contained in the basic equations is the second 
order one, two boundary conditions must be given for H, in order to determine 
the solution completely. Whereas he gave only one condition for H,, so that 
his condition is insufficient to solve the equation: generally. 

Moreover, there are some further inconsistencies in his paper. In the first 
place, he says that the total electric current in the direction Oy is zero. Ac- 
cording to his solution, however, it is not zero but finite, since 


total current ay : a2) dy = a {H,(co) — H 


fee. es ae 
Aa dy ,(0)} =e, : H,(0) si 0. 


0 


In the second place, he asserts that his velocity distribution reverts to the 
ordinary hydrodynamical one when S is zero. But this comparison is quite 
nonsense, because his solution is valid only when S is greater than unity. In 
fact, if we put S = 0 in his equation (27), the velocity distribution does not 
revert to the ordinary hydrodynamical one when Ry (= P,, in our notation) is. 
greater than unity. 

This paper has been written with an intention to comment on his paper and 
to propose an alternative solution. 


2. Fundamental Equations and Boundary Conditions 


The fundamental equations for the velocity of the fluid, v, and the intensity 
of the magnetic field, H, are derived from the modified NAVIER-STOKES equations 
of motion and MAXwWELL’s electromagnetic equations. Assuming, as usual, 
OxM’s law to be valid, we have the following dynamical and induction equa- 
tions for v and H in the steady state: 


ut 2 
(v- grad) v= — ry gradp + vy V2y + z : - (curl x Ay, (2.1) 
curl(y xX H)+7V7H=0, (2.2) 


together with the equations of continuity for the velocity and the magnetic 
fields: 


divy =0, (2.3) 


olive (0). (2.4) 


* 


ao 


<< h 
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) ie 


~~ 


where f is the pressure, o the density of the fluid, » and 7 (=1/4 0) the 
‘kinematic and magnetic coefficients of viscosity of the fluid respectively, and 
and o are the permeability and the electrical conductivity of the fluid res- 
-pectively. The electromagnetic units of system have been used throughout the 
present paper. 
E Now, we take the coordinate origin O at an arbitrary point on the wall, 
which is assumed to be electrically non-conducting and to have the same 
‘permeability as that of the fluid, and use Cartesian coordinates x, y and z with 
the axes Ox and Oy along and perpendicular to the wall respectively, and 
assume that 
| 1. all physical quantities are functions of the space coordinate y alone; 
_ 2. the velocity field is two-dimensional, that is, v= u(u,v,0) and the 
‘ suction is uniformly distributed over the whole surface of the wall; 
_ 3. the applied magnetic field is uniform and perpendicular to the wall, so 
that in the region considered, H = H(H,, H,, 0). 

The boundary conditions for the velocity field are exactly the same as 
those in ordinary hydrodynamics, namely 


p= O02 40, v= 0; yoo: u>U, (235) 


where U and v, denote the uniform x-wise velocity at far distance from the 
“wall and the suction velocity, respectively. 
5 Next, we shall consider the boundary conditions for the magnetic field. 
Since we have assumed that the wall is electrically non-conducting, no electric 
current flows in the wall. Therefore the magnetic field in the wall should be 
curl-free, that is curl H = 0. Thus, remembering the above assumption (1), we 
have —dH,/dy = 0 or H, = constant. This constant must be taken to be zero, 
since the applied magnetic field is normal to the wall. On the other hand, by 
‘virtue of the above assumption (1), divH = 0 implies that dH,/dy = 0 or 
H. = constant in the wall. This constant must be identical with the intensity 


y 
of the applied magnetic field. Therefore, we have 


4 H,=0, H,=H, (2.6) 


in the wall, where H, is the intensity of the applied magnetic field. 

: Now, both the normal and the tangential components of the magnetic field 
‘should be continuous at the fluid-wall interface, and the continuity condition 
‘for the normal component?) may be derived from div H = 0, while that for the 
| tangential component from curl H = 4 j by remembering that there flows no 


3) Strictly speaking, we should request the continuity of the normal component of the ‘magne- 
tic induction B’, instead of the magnetic field itself, since divH = 0 should be replaced by divB= 0. 
We have assumed, however, that the permeability of the wall is effectively the same as that of the 
fluid, that is, both the fluid and the wall are assumed to be non- magnetic. This assumption is not 
very restrictive and may be easily removed if necessary. 
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sheet current at the fluid-wall interface‘), where j denotes the electric current 
density. : 

On the other hand, at far distance from the wall, the magnetic field tends 
to the applied transverse field, since the x-wise velocity component ~ tends to 
the constant value U. Thus we have the following boundary conditions for the 
magnetic field: : 


y=0: H,=0, H,=H,; yoo: H,>0, H,>Hy. (2.7 


: 

It should be noted here that there is a constant electric field along the 
z-direction throughout the whole space. In fact, since curlH has the z-com- 
ponent alone, 7, = 7, = 0. Similarly, v x H has a component in the z-direction 
only. Hence, the combination of Oum’s law and the relation curlH = 4aj 
gives rise to E,=£,=0. Furthermore, since curl—=0 implies that 
dE ,/dy = 0, so that E, must be a constant, say Ey. Since the magnetic field at 
infinity tends to the contant value Ho, there is no electric current at infinity, 
so that the constant Ey must be equal to — w U Hy by virtue of OHM’s law. 
From the continuity condition of the tangential component of the electric 
field at the fluid-wall interface, we have the constant electric field —u U Ho 
throughout the whole space including both fluid and the wall. 

It should also be noted that the total flux of the electric current across the. 
direction Oy is zero, because 


lee) lo, 0} Le, °] 


io 1 t il GEL, 
total current =f 1, dy = rie / (curl H), dy = 9 | eres 
0 


v v 


0 0 
1 
This is one of the principal differences between GupTa’s result and the present 
one as has been mentioned in Introduction. 


3. Solutions and Shear-Stress 


Now we proceed to integrate the basic equations subject to the boundary 
conditions given above. 


By the assumption (1) in the preceding section, the continuity equation for 
the velocity field (2.3) becomes du/dy = 0 from which we have 


v= const = v%, (3.19 


where the constant has been determined by the boundary conditions (2.5) 
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Bestar, the continuity equation for the magnetic field (2.4) gives 
H, = const = Hy. (5.2) 


The fundamental equations (2.1) and (2.2) now reduce to 


0 dy —? aye + ang HO ay (3.3) 
Ce pe a 
ae * Sap Ta) =O (3.4) 
d @H, 
yt es) +n fhe (3:5) 


Re ke RL TT NS 6h eaceaamaae 


These three equations serve to determine three unknowns wu, H, and p. Induc- 
tion equation for H,, 1. e. 7 d*H,/dy? = 0 subject to the boundary conditions 
2.7), is auLomatically satisfied = virtue of (3.2). 


4 If we denote the’ pressure at infinity by ~,, it follows immediately from 
BS. 4) that 


j (oe Saee Ue go (3.6) 


“which shows the balance between the fluid pressure and the magnetic pressure 
_in the flow field, since H, is uniform throughout the whole flow field as shown 
"by (3.2). Since H,(0) = 0, the pressure at the wall is equal to p,. 

Integrating Equations (3.3) and (3.5) with respect to y and remembering 
that wu >U, dujdy >0; H,>0, dH,/dy > 0 at infinity, we have 


_, du ia err 
ade ~ reer Hy, H,+vU, (3.2) 
by ee Ao pees Ul He (3.8) 


dy 
Elimination of H, between (3.7) and (3.8) yields 


/ du Vo du v2 v gts 
| =f. PP ae eas Si 1 — 2 
a dy* » ul m) “a vn (1 S) u vn ( S) U , (3 9) 


“where P,, (= /y) is the magnetic Prandtl number and S (= ys Ho/4 2 @ vp) the 


'magnetic pressure number. Equation (3.9) is easily integrated, obtaining the 
following formal solution: 


=U Ae Be, (3.10) 


where A and B are constants of integration, and «, and a, are given as 


pee 1+ P,) + {2 ( +p — St a-5)}, (3.11) 


ey) 


~~ i 


+} 
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plus and minus signs corresponding to #, and «, respectively. It is found that — 
both «, and «, are always real, since 
2 2 

{2 1+ P,)} — 4° a — 5) =F {1 — Pt 44 SP} > 0. ; 

Let us first consider the case of uniform injection which corresponds to_ 
Vp > 0. In this case, «, and «, are both positive when S < 1, so that u diverges | 
at infinity unless both A and B are zero. Hence no steady solution can be : 
obtained in this case. On the other hand, when S > 1, «, is positive, so that the 
constant A should be zero. Then the remaining arbitrary constant is B alone — 
which would be determined by the boundary condition: (0) = 0°), so that there | 
remains no freedom to satisfy the boundary condition: H,(0) = 0. Thus we 
may conclude that in case when v, > 0, namely, the uniform injection is 
distributed upon the wall, no steady solution can be obtained for both cases of 
S <1 and S > 1 in the presence of “transverse magnetic field’. 

We next proceed to consider the case of uniform suction which corresponds 
to v, < 0. In this case, the expressions (3.11) for «, and «, may be rewritten as 
follows: 

a bee 


| 2 [a + Py) a5 ya 4 P,,)? - 4 P,, (1 Te s)] | 
Xs 


Il 


Vv 


= 5° [4+ P,) -VQ—P,2+4SP,]. | 


When S > 1, a, is positive, so that no steady solution can be obtained for the - 
same reason as in the case of vy > 0 and S > 1. Next, when S < 1, both of 
%, and %, are negative, so that the boundary condition at infinity for « can be 

satisfied for any values of A and B. Therefore, we shall deal exclusively with 

this case hereafter. 


Now the boundary condition at the wall for « requires that 
U+A+B=0. (3.13) 
Insertion of (3.10) into (3.7) leads to 
H, = Ff" [Wo— 9%) Ae + (09 — 9 4) Be™). (3.14) 
Since % and « are both negative, the boundary condition at infinity for Hy 


has been already satisfied. On the other hand, from the boundary condition 


») This case corresponds to Gurra’s solution. Thus, his magnetic field in the wall is not trans- 
verse but oblique to the wall. Moreover the angle between the magnetic field and the wall is a 
particular one determined by the value of H,,(0) (+ 0) given by his solution. 


gee ee ee 


: fol. XII, 1961 Hydromagnetic Flow over a Plane Wall with Uniform Suction 225 
at the wall, it is required that 

4 

4 (Vp —¥ a) A+ (v9 —va,.) B=0. ($515) 
Solving (3.13) and (3.15) for A and B, we have 

o 

A pe eames eet eee %. TF 

Z- v (%, — Og) 2a i 

3 (3.16) 
4 B de Sipehd 1 = Pat % oy | 

: Y (%, — 9) 2 o ‘ 


where « stands for /(1— P,)?+4S P,. 


: Insertion of A and B thus determined into (3.10) and (3.14) yields 

: uU 1 ey = 
4 iia =] aa Die ((1 Pi oo) 6 adie (1— Pi +0)2>"] , (3.17) 
and 

: Ay Pe U Pr Oy  ,ds¥ 

‘The above (3.17) and (3.18) together with (3.6) are the exact solutions for the 


hhydromagnetic flow over a plane wall with uniform suction in the presence of 
a transverse magnetic field, under the restriction of S < 1. It is to be noted 
that when S > 0 and/or P,, >0, u reverts to the ordinary hydrodynamical 
‘one obtained by GriIFFITH and MEREDITH [2]. : 

_ The related quantities such as the electric current density and the intensity 


“of the electric field may be obtained as follows. From the relation 


_ curlH 

| ot ee 

we have 

| es ie k Je p Renee 1 45 yy _ a pa keV) 
g.= 0, y= 9, ae a a) € (1+ P, + &) € mi) 

: (3.19) 
‘and again with the help of Oum’s law, we have 

~ oy ? 
.: == f] — = . 20 
: Boe Oty a+) Soa l, (3.20) 


“which has already been obtained in the preceding section. From (3.19) it will 
teadily be seen that 7, changes its sign at 


a In(j tt) 


me a Le Py 


: implying that the electric current flows in the negative z-direction in y < ¥, 
and in the positive z-direction in y > y,. The total current across the direction 
Oy is, of course, zero as was shown by (2.8). 


ZAMP XII/15 
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In the limiting case of very large P,,, the expressions for w, H, and7, become 
respectively | 


Uu 1 Bs “*.(1=S)y 2 
H, ye (3.21) 
~~  — é ‘ 7 
lal Uo 
a} 
; — Pe 
Ret me ie ~ O 
ty! deh ¥ ( : ), | 


excepting the very vicinity of the wall, implying that the electric current is 
confined in a narrow region near the wall and the induced magnetic field | 
changes steeply there. : 
Some typical distributions of u, H, and j. plotted against y for various 
numerical values of the parameters S and P,, are shown graphically in Figures 
1 and 2. 
The shear-stress at the wall, t, can be obtained from the formula: 


du 2 
t=or(5r) = me, (3.22) 


10 


Ue 


5-025 05075 ” 


Figure 1 


Figure 2 

Distributions of velocity, induced magnetic 

field and electric current density for various 
values of S when P,, = 1: 


Distributions of velocity, induced magnetic 
field and electric current density for various 
values of P,, when S = 0-5, 


ape Le al 
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a > 


which is identical with the one in ordinary hydrodynamics. Therefore we may 
assert that the shear-stress at the wall is not affected by the hydromagnetic 
interactions. On the other hand, the drag upon the wall due to the direct 
magnetic effect, i.e. the MAXWELL’s stress, is zero in the present problem, 
since (u Hy H,/4 2),-9 = 9. 


A The displacement thickness 6 may be obtained as follows: 
3 a u 
: b=| (ty) 4 
a 9 
4 “ (3723) 
, | Uo (1 ma S) (ee se 0) 
y 
bis Do ee, 0) , 


which shows that the displacement thickness increases with S but does not 
‘depend on P,, unless P,, is zero. It is quite surprising and somewhat paradoxical 
‘that the displacement thickness does not depend on P,,, since this means that 
‘no matter how small the electric conductivity of the fluid may be, 6 is affected 
by the presence of the magnetic field. Inspecting the limiting case when P,, is 
‘yanishingly small but not zero, we find that this paradoxical situation is due 
to the fact that an integral of the type 

ra 

fe e dy (e 20) 


0 
is not uniform with respect to the parameter e, namely, 


rd (4 (eee. 0), 
prety = 


0 


CO (e= 0). 


Then, in order to clarify the effect of the conductivity of the fluid upon the 
boundary layer thickness, we introduce tentatively an alternative thickness 
_6* defined by 


. 


co 


je Sf (1 is “Vy dy , (3.24) 


0 

which may be ‘a certain measure’ of the boundary layer thickness, although 
it has no direct physical meaning. 

By a simple calculation the value of 6* is obtained as 


? 1—S+ FP, 
ero Gas. Pay Pie 


t 
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Here — »/2 v9 (= 0%, say) corresponds to the ordinary hydrodynamical value, 
Thus, the ratio A of the thicknesses for the cases with and without the hydro- 
magnetic effect is given by | 


4q 


o* 1—S+ Ph 
= 5 6 } 
A= oe LR Wem EY (3.205) 
from which it is shown at once that 
One | Pin iy eee alae "= 
0S US SCE 9 ORG se 
aA Dal eres >0 A ees nes <0 4 : 
ost G—SPU+P,) =’ OP,  G-S\G+2,) = 


Therefore, we may say that the thickness of the boundary layer increases with 
P,, as well as with S, and the curves of A versus S and P,, would become 
concave and convex upward, respectively. 

It is to be emphasized here that the thickness 6* defined above may give 
only a certain measure of the boundary layer thickness. In fact, as may be seen 
from Figure 2 the curve of the velocity distribution for infinite P,, cuts those 
for finite F,, and approaches to unity more rapidly than those for finite P,,. 
Nevertheless, we may say at least that the boundary layer thickness in the 
case of non-zero P,, is thicker than that in the case when P,, = 0. 


m 


4. Energy Equation and Heat-Transfer 


In the steady state, the equation of energy may be written as: 
ec(v-grad) T—-AV?T = ee @, (4.1). 


where T is the temperature, ® the viscous dissipation, and c and 4 are the 
specific heat and the thermal conductivity of the fluid respectively. 
In the present problem, remembering the assumptions in Section 2, we have 


Ohi = Seger e © dH, du\? 
OQ c U9 dy A dy? _ => 16 we ( dy “) + O y a ) 4 (4. 2), 
Insertion of (3.17) and (3.18) into the above equation yields 
@T wv P. aT v8 P,U2 
dy= » dy B8cavi | 


» (4.3) 
* ne a ie es ac)? (1 = Fe =, a) es ae (1 ais Py at ax)? (1 7 ia + x) e* eet) | 


where P, (= ov c/A) is the Prandtl number. 
The solution of (4.3) subject to the boundary conditions: 


yes OF se 7 s yoo: Tl > T,, (4.4) 


SN 
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is given by 
S Too _ [1 OP te Py 2S Py gtePryl 
Ly — Too Hell t) Bt — a2 | 
, 2 Pee at ee ie ee | 
; OLS Soa BRR Eee ee 
P (+ Pn + %) (L = Pm +) etae)] 
B+o ‘ 


3 
where f stands for 1+ P,—,, and T,, and T,, denote the temperature at 
the wall and at far distance from the wall, respectively. 


The heat-transfer at the wall, g,,, can be obtained by the formula: 


| (Speco = 96% (Ty — Ts.) [1 — IE ae , (4.6) 


which is just identical with the ordinary hydrodynamical one. Thus, we may 
say that, just like the shear-stress, the heat-transfer at the wall is not affected 
by the hydromagnetic interactions. 


3 5. Concluding Discussions 

In this paper has been investigated the problem of the hydromagnetic flow 
over an infinite plane wall with uniform suction in the presence of a transverse 
‘magnetic field, and an exact solution for the steady state has been obtained 
‘when the magnetic pressure number is smaller than unity. It is concluded that 
both of the shear-stress and the heat-transfer at the wall are not affected by 
‘the hydromagnetic interactions but the boundary layer thickness increases due 
‘to the presence of the magnetic field. It has also been found that no steady 
‘solution can be obtained for the case of injection. 

Although we are concerned with the wall having infinitely large thickness 
in this paper, the problem may be easily extended to the case when a thin 
infinite flat plate is immersed edgewise in the flow field. In fact, assuming the 
flow field to be symmetric with respect to the plate, we would have the present 
solution for either half flow field in such a case. 

The solution obtained in this paper may be an asymptotic one in a sense 
that it may be realized on a long wall (or plate) with uniform suction at some 
‘downstream point from the leading edge. Evidently, the boundary layer begins 
‘to grow from zero thickness at the leading edge and tends asymptotically to 
the value given by (3.25), and the velocity profile and other quantities such as 
H, and 7, may attain the forms obtained in this paper. 

i The reason why no steady solution can be obtained in the case of S 2 1 
may be explained as follows. It is clear that the suction serves to check the 
growth of the boundary layer thickness so that the asymptotic steady state can 
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be realized. On the other hand, the presence of the magnetic field counteracts 
this checking effect of the suction, so that the asymptotic thickness of the 
boundary layer increases with S and ultimately when S (= we H3/4 2 0 vo) 
reaches to unity the counter effect of the magnetic field balances with the 
checking effect of the suction. Therefore no steady solution can be obtained for 
S = 1, just like in ordinary hydrodynamics no steady solution exists in the 
case of no suction. According to this argument it is quite natural that no 
steady solution is obtained for the case of S 21, and also for the case of 
injection. 
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Zusammenfassung 


Diese Abhandlung behandelt die Str6mung einer inkompressiblen, zahen und 
elektrisch leitenden Flissigkeit tiber eine ebene Wand mit homogener Absaugung. 
und senkrechtem Magnetfeld. Eine exakte Lésung der hydromagnetischen Glei- 
chungen kann erhalten werden, wenn die magnetische Druckzahl kleiner als eins 
ist. 

Aus den Rechnungen findet man, dass sowohl die Reibungsschubspannung als 
auch die Warmeiibertragung an der Wand unabhingig sind von der hydromagneti- 
schen Wirkung, aber dass die Grenzschichtdicke zunimmt mit der Zunahme des. 
hydromagnetischen Einflusses. 

Es sei bemerkt, dass keine stationdre Lésung erhalten wird fiir den Fall des 
Ausblasens. 


(Received: September 20, 1960.) 
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‘ Rotationssymmetrische, wirbelbehaftete Stromung 
als Variationsproblem 


s Von BoHDAN KRAJEWSKI, Warschau, Polen 


Einleitung 


Die stets weitergreifende Verwendung der Str6mungsmaschinen bedarf 
einer eingehenden Kenntnis des Geschwindigkeitsfeldes, die es erlaubt, hieraus 
die richtige Schaufelform zu ermitteln. Die bekanntesten Arbeiten auf dem 
Gebiet der Berechnung der Strémungsmaschinen sind die von RuDEN [1}}), 
/ Wu [3], GOLDSTEIN [4], MONROE [5], STRSCHELETZKY [7] und GRAVALOs [8]. 
‘In den bis jetzt veréffentlichten Arbeiten verwendet man vor allem die Itera- 
-tionsverfahren zur Lésung auftretender Differentialgleichungen. In der vor- 
_liegenden Arbeit wird ein neues Verfahren zur rechnerischen Ermittlung des 
_Geschwindigkeitsfeldes vorgeschlagen, wobei die Begrenzungen der Strémung 
sowohl aus festen Wanden als auch aus Stromflachen bestehen kénnen. Die 
Form dieser Begrenzungen entspricht dem Zustand des gasdynamischen 
-Gleichgewichtes des strémenden Gases in dem betreffenden Raum. Die An- 
-nahme der rotationssymmetrischen Str6mung erméglicht zum Beispiel — zwar 
-angendhert, aber meist mit geniigender Genauigkeit — die in Turbomaschinen 

vorkommende Strémung zu berechnen. In dieser Abhandlung wird sowohl die 
absolute Bewegung, die im Leitrad und im Ubergangsraum vorkommt, unter- 
‘sucht als auch die relative Bewegung, die im Laufrad stattfindet. Im folgenden 
werden also die wirbelbehafteten Str6mungen untersucht, deren Geschwindig- 
-keitsfelder axiale, radiale und Umfangskomponenten der 6rtlichen Geschwin- 
-digkeit aufweisen, wobei angenommen wird, dass sowohl die Verteilung der 

értlichen Geschwindigkeiten, Driicke und Temperaturen als auch die Verteilung 
der totalen Enthalpie nur von der axialen und radialen Koordinate abhangen. 


1. Die Grundgleichungen 


Fiir die stationare Strémung eines reibungsfreien idealen Gases konstanter 
“spezifischer Warme gelten bekanntlich die Gleichungen: 
-Zustandsgleichung 

Va hT; (1) 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 241. 


DV BouDAN KRAJEWSKI : ZAMP 
Kontinuitatssatz dieoo 0: (2) 
Impulssatz 
De af Ge 1 : 

= Of eS 2 pe d : 3) 

= = grac($) c x rote » era p (3) 

Energiesatz : 
T (>) = €erad S = 0), (4) 


wobei vorausgesetzt ist, dass in der Strémung keine Energiedissipation und 
keine Warmeleitung stattfindet. Die Gleichung besagt dann, dass die Entropie S 
in Richtung der Stromlinie bzw. auf der ganzen Stromflache einen konstanten: 
Wert hat. Das bedeutet, dass die Entropie eine Funktion der Stromfunktion 
sein muss: 


S=f(%), (5) 
Nach dem ersten Hauptsatz der Warmelehre gilt bekanntlich 


dh= Tas +? (6) 


Durch Einfiihren der Gleichung (4) in die letzte Beziehung erhalt man 


DW © A ADE : 
(Br) - +) 7 

Fiir die totale Enthalpie gilt folgende Abhangigkeit : 
H=h+. (3) 


Bei Beachtung der Gleichungen (3) und (7) erhalten wir fiir die Energiegleichung 


den Ausdruck 
DH 
“pyr Ole rerot chaz. (9) 


Dabei ergibt sich wegen ahnlicher Zusammenhange wie bei der Entropie S 
Ve ee a (10) 


Nach der Elimination von h und # aus den Gleichungen (3), (6) und bei Beach- 
tung von Gleichung (8) erhalten wir als grundlegende Beziehung fiir die weitere 
Untersuchung 

grad H — T grad S = c x rotc. (11) 


Diese Gleichung, die zuerst von VazsonyI [2] hergeleitet wurde, spricht den 
verallgemeinerten Croccoschen Wirbelsatz aus und bezieht sich nur auf die 
absolute Strémung, die im Leitrad bzw. im Ubergangsraum vorkommt. Um die 


aT 


% 
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entsprechende Beziehung fiir die relative Strémung herzuleiten, fiihren wir den 
Begriff der relativen, totalen Enthalpie ein. Es ist nadmlich 


2 2 
Jah+ 5-4 =H-o(C). (12) 


Auf Grund der Beziehung fiir die absolute Geschwindigkeit 

} c=wu, (13) 
wobei 
, u=OxXxr ist, (14) 


erhalten wir fiir den Impulssatz die Gleichung: 


De _ Dw 1 
P= pr tt Zax w=—~— gradp, (15) 
wo 
Dw w? 
| a sae grad (+) — wx rotw (16) 
bedeutet. 


Bei Beriicksichtigung der Gleichungen (6), (12), (15), (16) und nach einer 
kleinen Umformung ergibt sich als grundlegende Beziehung fiir die relative 
Str6mung 

grad J — TgradS= wx rotw+wx2@. (17) 


Verwendet man die Gleichungen (13) und (14), so erhalt man 


rotc =rotw+ 2a. (18) 


Nach Einsetzen von (18) geht die Gleichung (17) tiber in 


grad J — T grad S = w x rotc. (19) 


Dieser Zusammenhang ist zuerst von Wu [3] angegeben worden. 


2. Krummlinige Orthogonalkoordinaten 


Um die vektoriellen Gleichungen (11) und (19) in drei skalare Gleichungen 
zu zerlegen, fiihren wir ein lokales Bezugssystem ein. Wahlen wir als Parameter 
den Drehwinkel @ in der Meridianebene, in welcher die Parameterlinien ortho- 
gonale Netze bilden, dann ergibt sich nach Drehung um die Drehachse x die 
Darstellung der Rotationsflachen. Wie zuvor erwahnt wurde, hangt das Ge- 
schwindigkeitsfeld der Strémung nur von der axialen und radialen Koordinate 
ab, und deswegen geniigt es, die Strémung in der Meridianebene zu betrachten. 
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Fiir unseren Fall, wo es sich um orthogonale Koordinaten handelt, wird zweck- 
miassig der Winkel  eingefiihrt, der im allgemeinen eine Funktion des Ortes ist 
und von x, y oder m,  abhangt (Figur 1). Werden die értlichen Kriimmungs- 


x 


Now 2) 
Figur 1 
Krummlinige Orthogonalkoordinaten. 


radien der Meridianlinie und der Normallinie entsprechend mit o,, und , 
bezeichnet, dann lassen sich die zugehérigen Verschiebungen dm und dn langs 
der Meridianlinie und Normallinie angeben durch 


dm = Om AP , dn rs On ap, : (2.1) 


Im folgenden wird die entsprechende Kriimmung in M als stetig und von null 
verschieden vorausgesetzt. Fiir ein Bogenelement gilt bekanntlich 


ds =dmi,,+rdOig + dni, , (2.2) 


i, 1, bedeuten die Einheitsvektoren, die entsprechend in die Richtung der 
Tangente der Kurve «m» oder «n» fallen. ig bedeutet den Einheitsvektor in der 
Umfangsrichtung. 

Die Vektoren i,,, ig sollen mit i, ein Rechtssystem bilden. Nach Einsetzen 
von (2.1) in die Gleichung (2.2) geht diese iiber in 


ds = Om IP i, 1 Y dO ig ate On dp, ie - (23 3) 


Im Falle krummliniger orthogonaler Koordinaten lasst sich der Wirbelvektor C 
in folgender Weise als Determinante aufschreiben 2) 


| Om tn , Y ig ’ On i, , 
1 0 7) 0 
ae eee ew ie (2.4) 


m7 On | 2m’ 08° dp,’ 
Om Ce ? fe (es ? On C . | 


®) Vergleiche [6]. 
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Bezieht man die Analyse der Strémung im betreffenden Falle auf die ortho- 
‘gonalen Koordinaten m, n, so ergeben sich Zusammenhiange, welche die weitere 
Untersuchung wesentlich vereinfachen. 

- Da die Strémung rotationssymmetrisch ist, verschwinden die Ableitungen 
nach @. Ausserdem ist in der Meridianebene m, » nur die Meridiankomponente 
der Geschwindigkeit C,, vorhanden, wobei C, = 0 ist. Unter Beriicksichtigung 
‘des vorhergesagten und nach Einfiihrung der Operatoren 


2 i) 0 0 i 
5 On 9Pn OfmeOs OO om ’ 


ergibt sich anstelle der Gleichung (2.4): 


Se PPO Gye Om Cm) = O(r Cy) 
a yon 'm” Om ON free Yer CS 
Verwendet man die Beziehungen 
OH OH. 0H 
| ie 088" On.” | 
und (2.6) 
; 0S 0s Os 
ere? im ™! x0@ 8." on,” | 


so ergeben sich nach Zerlegung der Gleichung (1.1) in die drei Richtungen 
4, ig, 1, die drei skalaren Gleichungen 


dH 7 0S _ Cy Or Cy) 


Om om Y om , 2) 

OH OS ig OEY) 
SOGt eae iawn.” ae) 
eet os ile, Cala lune Gd ae 
“On on Cm Om On ser on * ee 


Aus fritheren Betrachtungen (GI.(4) und (9)) geht hervor, dass sowohl die 
Entropie S als auch die totale Enthalpie H auf ein und derselben Stromflache 
einen konstanten Wert hat. Da jeder Meridianlinie nur eine Stromflache zuge- 
ordnet ist, ergibt sich dann aus der ersten Gleichung: 


oH Fae eee 00 'Cy). -. 
om om 


Fiir die zweite Gleichung erhalt man wegen der Rotationssymmetrie die Be- 
ziehungen: 
OH OS) 


mee se 
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Daraus folgt, dass die Gleichung (2.9) die Grundgleichung der rotationssym- 
metrischen, wirbelbehafteten Str6mung darstellt. Zur weiteren Umformung der 
letzten Gleichung gehen wir aus von der Darstellung des Wirbelvektors in 
Zylinderkoordinaten. In diesem Falle lasst sich der Wirbelvektor folgender- 
massen als Determinante aufschreiben: 


i j k | 

if 0 i) 0 4 

6 =rote= ian? nee ae (2. 10) 
Ge nie, c.. | 


Da nun die Einheitsvektoren ig und j sowohl zu der m, n-Ebene als auch zu der 
x, y-Ebene senkrecht gerichtet sind, miissen die Umfangskomponenten des 
Wirbelvektors gleich sein: 

DeCa sbee, acy 


Om On or vas Co | 
Setzt man die letzte Beziehung in Gleichung (2.9) ein, so gilt 
OH Os dC. C, 9(r C,) 
“On Lae On Ce (% Or he ol ar Bie eS (21 


Bei der rotationssymmetrischen Strémung ist es méglich, eine Stromfunktion 
zu definieren. Durch Einfiihrung einer Stromfunktion lasst sich die Kontinui-. 
tatsbedingung identisch befriedigen. An die Stelle der Stromlinien treten dort 
die Stromflachen, die Rotationsflachen sind. Jeder Stromflache kann ein den 
Durchsatz charakterisierender Parameter Y zugeordnet werden. Der Betrag_ 
des Gradienten ist bei rotationssymmetrischer Strémung gleich der mit 7 multi-. 
plizierten Stromdichte. Es ist namlich | 


o7C,, = 09 | grad | . (2.13) 
Daraus ergibt sich fiir unseren Fall 
Ones 
or C. == '06 7% : (2.14) 


Da die Meridianlinien «my» eine einparametrige Kurvenschar derart bilden, 
dass durch jeden Punkt der Meridianebene eine m-Kurve der Schar hindurcht 
geht, ist es nun méglich, folgenden Operator einzufiihren: 


0 0 
0) 0 er, on’ (2:15) 


Setzt man die eben erhaltene Beziehung in Gleichung (2.12) ein, so erhalt 
man nach einer kleinen Umformung 


re OH rep AS 9 MWC) 2) , WC, 06, 
0 OV 20 ov a Ow oy (eer (2.16) 


: aw aw 

¥ e7C,= >,» e7C,=— 5 >> (2.17) 
so geht die Gleichung (2.16) iiber in 

: 

74 re 0H ve 79S _ _@ ArC,)*/2] | A Qo (Ov 

: Oo OW o Ow Y Qo Ow or Lre ( or )| (2.18) 


0 ow 

i ee er a Ox )]- 
Tn Analogie zum letzten Vorgehen kénnen wir die entsprechenden Beziehungen 
auch fiir die relative Bewegung herleiten. Es wird namlich 


5 Pipe an ES ae Wy or Cy) 
4 On on Or Ox Y On 


(2.19) 


Es kommen noch zusiatzlich die Beziehungen hinzu, die die Gleichungen (13) 
und (14) liefern, namlich (Figur 2): 


ren oO = to ae 


u 


Figur 2 


Das Geschwindigkeitsdreieck. 


Nach Einsetzen der letzten Beziehungen und bei Beriicksichtigung von Glei- 
chung (2.15) ergibt sich dann nach kurzer Umformung 


re OJ ro 7S __ reo HC) , @ Aly Cy)*/2} 
Bae. 00 ov 00 OV Y Oo Oy 


ter ieete lt oleate) 22) 


3. Variationsproblem 


Wir gehen jetzt zur Untersuchung des Variationsproblems fiir ein Doppel- 
integral iiber: 
i ay) ov 


BLP = [lil res, #) : ica | dx dy = Extremum. (3.1) 
B 


| 
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Es wird die Funktion ¥(x, 7) gesucht, die das Integral F zu einem Extremum 
macht und am Ende des Bereiches « B» gegebene Werte annimmt. Das Integral 
F besteht nur dann, wenn die Eulersche Differentialgleichung des Variations- 
problems erfiillt ist. Fiir unseren Fall lautet diese Gleichung 


Ae ee 3 laye) =o (3.2) 


Um unterscheiden zu kénnen, ob eine Extremale wirklich ein Extremum 
liefert, muss zusadtzlich die Legendresche Bedingung erfiillt werden. In unserem 

Falle lautet diese Bedingung | 
of 

OV, 


cecbl : 
+03 Fy abogpiek ee | 
Es ist leicht zu priifen, dass die Gleichungen (2.18) und (2.20) die Eulersche 
Bedingung (3.2) erfiillen, so dass die zugehGdrigen Integrale sich folgender- 
massen schreiben lassen: 


1. fiir die absolute Bewegung 


F@ “ale (ie [z() — TSH) - Ss) 


+ Sel (RY + (SE) T} eee evomm, 68 


mit lar C5 


2. fiir die relative Bewegung 


Fey = ff {72 [7 — 7 SH) +0 re) — 0) 


stig oO (S-) aie (s-) |} dx dy = Extremum. (3.5) 


Die letzten beiden Gleichungen sind auch einer physikalischen Deutung fahig. 
Wir betrachten zu diesem Zweck das Integral (3.4). Bei Beriicksichtigung der 
Beziehungen (6) und (8) erhalten wir: 


dT TUS a(< | Bae 
0 
wobei 
C2 Ca C2. 


Verwendet man zusatzlich die Beziehung 


Gi Oa BCE Re) 


ra 
f 
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dann ergibt sich anstelle der Gleichung (3.4) 


: FO lla \r 0 [ez a \ dx dy = Extremum . (3. 6) 


Bezeichnen wir die Masse @ C,,, welche durch die Flacheneinheit eines Quer- 
‘schnittes in der Zeiteinheit fliesst, als die Stromdichte, dann kann die letzte 
Gleichung wie folgt gelesen werden: Liegt die rotationssymmetrische Str6mung 
fest, dann erreicht die Summe der Impulsdichte 


4 d 
ces glia 
C(¥) hat einen konstanten Wert auf der ganzen Stromflache] 


innerhalb des Rotationshohlraumes einen Extremwert (Minimum). 


4. Randbedingungen 


Aus den vorangehenden Betrachtungen geht hervor, dass jeder Strom- 
flache ein den Durchsatz charakterisierender Parameter WY zugeordnet ist. 
Dieser Parameter andert seinen Wert beim Ubergang von einer Stromflache zu 
der anderen, stellt aber eine stetige Funktion sowohl innerhalb des Raumes als 
‘auch an seinen Grenzen dar. Wahlen wir nun die Kontrollflache, die senkrecht 
zu der Symmetrieachse steht, dann geht sie durch jede angetroffene Strom- 
flache hindurch. 

- Die stetige Funktion stellt nun die Durchsatzverteilung langs der Linie 7 
dar und nimmt am Rande des Integrationsbereiches die vorgeschriebenen 
Werte an (Figur 3). Denken wir uns die gesuchte Funktion (x, 7) durch den 
folgenden Ansatz angenahert: 


W(x, 1) = (x, 7) + O(%, 1, 41, Ay, «++ » Ay) (4.1) 


Kontrollflache 


> 

x 
Figur 3 

Die beiden Koordinatensysteme (x, y) und (m, n). 
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Dabei ist p(x, 7) eine beliebige Funktion, welche die angenommene Durchsatz- 


verteilung darstellt und folgenden Randbedingungen geniigt : | 
G | 

Pol%,%:)=9; Pol%, 72) = Ve = W(x,7,) , (4. 2) 

wo | 


G[*] verfiigbarer Durchsatz; 
r(x) Radius der inneren Begrenzungswand; 


7 (x) Radius der dusseren Begrenzungswand. 


Dagegen lauten die zugehérigen Randbedingungen vom zweiten Glied des An- 
satzes (4.1): 
PUG, His bys Gano 2 ip hy) = O(%, Ves Cys Wy Seay a) Oy (4.3) 


wobeil 4, a, ...a, noch zu bestimmende Konstanten sind. Benutzt man die: 
Bezeichnungen gemass Figur 3, dann gilt fiir die Funktion g(x, 7): 


G  — #2(x) 
ne) = ee (4.4) 


Yo r(x) = v; (x) 


Dagegen asst sich die Funktion (x, 7, a, a, ...a,), die die gesuchte Durch- 
satzverteilung darstellt, folgendermassen schreiben: 


Hx,4, @y, Ay, -..,@,) = (7,(%) =r] [r — r;(x)] (@, + 4px + agr+ G@gX¥...). Taam 


Wenn die a; so bestimmt werden sollen, dass F(W) ein Extremum wird, so. 
miissen die partiellen Ableitungen von F nach a, verschwinden 

ai ger, ES eats 
0a, 0a, 0a, 


= 0, (4.6) 


Mit den aus diesen Gleichungen sich ergebenden Werten der a; gehen wir in (4.1) 
ein, womit die gesuchte Naherungsfunktion gewonnen ist. 


5. Rechnerisches Vorgehen 


Um eine numerische Lésung von Gleichung (3.4) und (3.5) gewinnen zu 
k6nnen, miissen noch weitere Annahmen getroffen werden, die die rechnerische 
Ermittlung wesentlich vereinfachen. Zu diesem Zweck wird unsere Betrachtung 
nur auf solche Funktionen H(¥W), S(W), [(W) beschrankt, die linear von YW 
abhangen. Es kénnen natiirlich Bedenken entstehen, ob nicht viele technisch 
wichtige Probleme bei einer derartigen Beschrankung ausser acht gelassen 
werden. Jedoch lasst sich die Mehrheit der in der Technik auftretenden Pro- 
bleme erst mit dieser Annahme behandeln. Auf Grund dieser Voraussetzung 


, 
: 


« 
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; liefert uns die Gleichung (4.6) «7»-lineare Gleichungen mit den «2» Unbekann- 
_ ten a, dg, 3, ... , 4,, die sich in tiblicher Weise errechnen lassen. Es steht natiir- 
; lich nichts im Wege, von der Potenzreihenentwicklung der Funktionen H(¥), 
_ S(W), /(W) noch mehr Glieder in Betracht zu ziehen, aber die rechnerische 
_ Ermittlung wird dadurch erheblich erschwert. Demgemass ergeben sich fol- 
gende Beziehungen: 


Heakiphgl, S=m,4+ m2, T= n+ te, (5.1) 


Mero ee 


wo k, m,n bekannte, sonst beliebige Konstanten sind. Die Wahl dieser Kon- 
stanten ist von der gestellten Aufgabe abhangig, es kann aber frei tiber sie 
verfiigt werden. In unseren Integralen (3.4) und (3.5) treten noch die Funktio- 
nen g,/0 und g/o, auf, die am Anfang der Rechnung unbekannt sind. Um diese 
_ Schwierigkeit zu beseitigen, benutzt man zuerst die Werte von @,/e und @/G, 
- die sich aus der eindimensionalen Betrachtung ergeben. Die Voraussetzungen 
- seien die gleichen wie beim Vorentwurf, das heisst, wir ersetzen die komplizierte 
_ Geschwindigkeitsverteilung in den betrachteten Kontrollquerschnitten durch 
- Mittelwerte. Mit den auf diese Weise gewonnenen Werten bilden wir entspre- 
- chende Ersatzfunktionen von @,/o bzw. g/g) und gehen dann in (3.4) baw. (3.5) 
ein, womit die gesuchte Naherungsfunktion gefunden ist. Man konnte die 
_ Rechnung noch verfeinern, indem man mittels der gefundenen Funktion Y 
eine neue Ersatzfunktion fiir 9/0) und ,/@ bildet und die ganze Rechnung 
wiederholt. Der erste Rechnungsgang liefert jedoch schon eine geniigende 
Genauigkeit fiir die in der Technik vorkommenden Probleme. Die Geschwin- 
digkeitskomponenten lassen sich aus folgenden Beziehungen berechnen : 


ep ee Gras 2 Orn ee La aS (5.2) 


a @ ot. o fox” 2 Y 


Hat man das Geschwindigkeitsfeld ermittelt, dann sind wir imstande, die 
Schaufelform zu bestimmen. 
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Summary 


(amongst other things) in turbomachinery. The problem can be reduced to a partial 
differential equation with given boundary conditions. Bearing in mind that the 
differential equation considered corresponds to a variational problem, the Ritz | 
method is used for a numerical solution. This requires a few additional assumptions _ 
to facilitate the computation procedure and enables us to correlate various design 
parameters. The velocity field being thus obtained, it is possible to determine the 
form of the blades. 


; 

3 

This is an analysis of the axially symmetric rotational flow which takes place ; 
: 


: 


(Eingegangen: 9. Juni 1960.) 


Treftertheoretische Untersuchung von radioaktiven 
Einwirkungen auf geléste Substanzen’) 


Von Kari-HEtnz Miter, Saarbriicken?2) 


Einleitung 


Durch Bestrahlung treten oft bei biologischen Objekten entscheidende 
Umwandlungen auf, die auch auf benachbarte, nicht unmittelbar der Strahlung 
ausgesetzte Bereiche des Organismus iibergreifen. Man spricht in Anbetracht 
dessen von einer primaren Wirkung, der durch die Strahlung an empfindlichen 
Molekiilen statthabenden Umwandlung, und von einer sekundiren Wirkung, 
der Summe aller Effekte, die eine primare Umwandlung mit sich bringt. Die 
wohl wichtigste Sekundarreaktion ist die Diffusion. Durch sie kénnen eine 
Reihe weiterer in Gang gebracht werden. Man denke dabei an das Heranschaffen 
von Schutzstoffen gegeniiber Strahlenwirkungen — ein Positivum —, daneben 


auch das Infizieren «gesunder» Zellen durch umgewandelte Molekiile — ein 
Negativum. 


1) Geférdert durch das Bundesministerium fiir Atomenergie und Wasserwirtschaft, Bonn, 
Kontrakt: FN 1.222, 


*) Universitat des Saarlandes, Institut fiir Mechanik. 
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das Erfassen der sich im einzelnen abspielenden molekularen Prozesse (Treffer) 
zu analysieren. Unter einem Treffer versteht man hierbei die Aufnahme einer 
fur die jeweilige Molekiilart charakteristischen Strahlungsdosis, die eine Um- 
_wandlung bewirkt. 
Wahrend die «klassische» Treffertheorie nur mit relativ speziellen Reak- 
‘tionen arbeitet, namlich solchen, bei denen das Wirkvolumen von der bereits 
‘aufgenommenen Dosis unabhangig ist, und mégliche riicklaufige Prozesse un- 
beachtet lasst, soll bei uns auch auf Vorgange dieser Art «treffertheoretisch» 
-eingegangen werden. Wir wollen jedoch die oft beobachtbare Tatsache, dass 
_ diffundierende Teilchen Veranderungen an anderen Teilchen verursachen, nicht 
periicksichtigen. Sekundarreaktionen solcher Art lassen wir ausser acht; wir 
‘untersuchen also nur Direktumwandlungen bei unmittelbar auftreffender 
- Strahlung und spontane Umwandlungen (zum Beispiel Erholung). Die Abhan- 
gigkeit des Wirkvolumens von der bereits absorbierten Strahlungsdosis vermag 
-man auch zu erfassen, wenn zwischen dem «Ausgangszustand» («) — Molekiil 
vor der Bestrahlung — und dem «Endzustand» (w) — Molekiil nach vollzogenem 
_ Treffer g— ein oder mehrere «Zwischenzustande» (B) (y) ... eingeschaltet werden, 
das heisst solche, bei denen nur ein Teilquant (gp, 7,, ...) der zum Erreichen des 
Endzustandes erforderlichen Strahlendosis appliziert ist. Schematisch sieht 
- dies so aus: 
; (x) + 9, > (A), | 
(8) + 4, >)» 


S (o —1) +4, > (@), 
_ wobei sich insgesamt die Reaktion abspielt: 


(x) +9 > (@); a= 30, (2) 


Gegenlaufige Spontanreaktionen konnen in den einzelnen Stufen des Schemas 
_ (1) auftreten, wenn die Lebensdauer kiirzer ist als der Zeitraum bis zur Auf- 
- nahme des nachsten Quants. 

Ein Teil der einfallenden Strahlung wird von der nicht-strahlungsempfind- 
lichen Substanz inklusive den Molekiilen des Endzustandes absorbiert und so 
der radiochemischen Wirkung entzogen; dies bedeutet oft einen gewissen 
- Schutz fiir die Zellsubstanz selbst. 

Wir wollen unsere Untersuchungen in folgende Abteilungen gliedern: 

1. Stérungsprozesse, das heisst Prozesse wahrend der Bestrahlung, 

2. Ausgleichsprozesse, das heisst Prozesse nach der Bestrahlung. 


Wahrend die Ausgleichsprozesse allein durch die Diffusionsgleichung regiert 
- werden, benétigt die Beschreibung der Stérungsprozesse nichtlineare partielle 


| 
| 
: 
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Differentialgleichungssysteme, deren Integration in der Regel grossen mathe- 
matischen Aufwand verlangt. Um diesen Aufwand nicht ins Uferlose wachsen > 
zu lassen, beschranken wir unsere Betrachtungen auf relativ einfachgebaute : 
Strahlenquellen wie | 
a) den Parallelstrahler (Parallelstrahlen iiberall gleicher Flachendichte: Ebene) ; » 
b) denachsensymmetrischen Strahler (Strahler ist geradlinig: Draht, Réhrchen) 
und | 

c) den punktsymmetrischen Strahler (Strahler ist punkt- oder kugelférmig: 

Kugel). : 

Diese Typen erfassen die in der Physik und Strahlentherapie am haufigsten — 
auftretenden Falle, stellen also fiir die Praxis keine allzu fremden Spezifika- 
tionen dar. 

Der Vereinfachung halber soll die Reibung des Zellinhalts an den Zell- 
wanden nicht in die Rechnung eingehen, auch nicht Temperaturschwankungen | 
und Inhomogenitaten des Lésungsmittels, in das die strahlungsempfindlichen 
Molekiile eingebettet sind; fernerhin soll die Wahrscheinlichkeit, von einem 
Quant getroffen zu werden, fiir alle Molekiile des gleichen Zustandes gleich 
sein. Wir arbeiten also mit Mittelwerten. Die zugehérigen Diffusionsgeschwindig- 
keiten seien konstant. 


1. Sté6rungsprozesse 


Der ablaufende Prozess setze sich aus einer Reihe von Teilreaktionen zu- 
sammen. Durch die Absorption eines Quants werde ein Teilchen des Ausgangs-. 
zustandes («) in den Zwischenzustand (8) und durch Absorption weiterer 
Quanten Schritt fiir Schritt in die iibrigen (m — 1) Zwischenzusténde umge- 
wandelt. Nach Aufnahme eines Quants durch ein Molekiil im letzten Zwischen-_ 
zustand werde der Endzustand (w) erreicht. Solchen Direktumwandlungen 
k6nnen nun auch gegenlaufige Spontanreaktionen tiberlagert sein. Beriicksich- 
tigen wir auch solche Vorgange, dann gehen neben den Wirkungsvolumina 
(vq, ¥ -.. ¥) und den Diffusionskonstanten (%2,%p, ++. ,,) noch die sogenannten 
Riickumwandlungskonstanten (Up, My, --. 4) als weitere Parameter in die Dif- 
ferentialgleichungen ein. Wir bezeichnen fiirderhin die Flachendichte der Strah- 
lung mit J’ = J(t, t); r ist dabei der Ortsvektor. Auf Grund der Poisson-Stati- 
stik, deren Giiltigkeit wir hier voraussetzen wollen, und den Gesetzen der Dif-_ 
fusion errechnet sich die zeitliche Anderung der Konzentration N i(t, 4), das 
heisst der Anzahl der Molekiile pro Kubikzentimeter des i-ten Zwischenzu- 
standes, zu: 

ON; 
aie #; AN; + (v1 Nia — 9; Ni) J+ isa Nia — oN) | (3) 
(l<i<m), | 


(fiir Ausgangs- und Endzustand gelten Differentialgleichungen gleicher Art). 


oe 
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_ Das so gewonnene Differentialgleichungssystem ist wegen der Produkt- 
pildung der zunadchst unbekannten Konzentration N mit der Flachendichte J 
nicht linear. Auch ist das System nicht vollstandig, es fehlt noch eine Aus- 
‘sage iiber die Variation der Strahlungsdichte J. Die Tatsache, dass der fiir den 
Quantenstrom undurchlassige Teil des Raumes proportional der absorbierten 
‘Strahlung ist, fiihrt uns zur fehlenden Relation zwischen J und den méglichen 
_(m + 2)-Zustanden N;: 


4 OF s. 
; Or a J 


(4) 


A m +2 
: + 2% N | ss 
Es ist dabei |r| =7 gesetzt. A misst die Divergenz der Strahlung und die 
Absorption durch das Lésungsmittel (A = Absorptionskoeffizient [cm™)) : 


Strahler Q | A 

Parallelstrahler ue (5) 
Achsensymmetrischer Strahler 1+Ar | 
Punktsymmetrischer Strahler 2+Ar 


Zu Beginn der Bestrahlung sei: 


| No= const 
i= 0: Bon ho 


6) 
fen 


Das Verhalten der Molekiile an der Oberflache r = ry des bestrahlten Kérpers K 
(Figur 1) erfassen wir mit: 


oN; 
Ee h.N, |. =F, 7 
fa, Gt + 8Ni] =F (to), (7) 

F; 
a To 

K 
a 
fp 
Figur 1 


Ein Teil der Strahlung trifft auf einen strahlungsempfindlichen Kérper. 


_ wobei 0/0 die Normalableitung, a; und 6; Konstanten und F;(ro, t) vorgegebene 
- Funktionen bezeichnen. 


> 


~~ ] 
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Ein Mass fiir die den Kérper K zwischen P, und P, treffende, also «ein- 
fallende» Strahlung, ist: | 
La hk Ake (8) 
Diffundieren nun noch Schutzstoffe von aussen her in K hinein, dann ist deren. 
Konzentration ebenfalls durch (3), (6) und (7) erfassbar. Zu beachten ist dabei 
lediglich das Verhalten der Molekiile dieser Stoffe im Innern von K, also die 
Kenntnis der zugehorigen Parameter x, A, yu, v. | 
Wird jedoch die Anzahl der reagierenden Molekiile im Kérper (vom Volu- 
men V) nicht gedndert, dadurch etwa, dass die Wandung fiir sie impermeabel 
ist oder dass die gleiche Anzahl ein- wie austritt, dann gilt noch 


, 
: 


[Ena anv. (9) 


t=1 


Die alle Vorgange beschreibenden Gleichungen (3-9) erfahren nur, das sei noch 
kurz erwahnt, eine geringfiigige Modifikation, wenn man statt mit den Konzen- 
trationen N; mit den Wahrscheinlichkeiten W, fiir das Zustandekommen von 
Trefferereignissen rechnet ; wir brauchen deshalb hier nicht naher darauf einzu- 
gehen. 

Zum Schluss der grundsatzlichen Betrachtungen wollen wir noch die Spalt- 
prozesse von Molekiilen erwaéhnen. Die bisher angegebenen Gleichungen sind 
hierbei nun auch auf die Spaltprodukte anzuwenden, und zwar auf folgende 
Weise. Zur Hlustration wahlen wir ein einfaches Beispiel, den «Eintreffervor- 
gang» mit einfacher Spaltung. Ist N, die Konzentration des Ausgangszustandes, 


Ng und N, die der Bruchstiicke, dann gilt: 


ON ONp oath 
Bae AN ch PING] Sie RRS Pete es 
ON, 


- 1 oy A 
HE MAN, — az %—eNat, Ge =I [A +y,NetogNg tr, N,]. 


Vv 


Synthetisierungsvorginge erfasst man auf ahnliche Weise. 


2. Ausgleichsprozesse 


Auch eine grosse Anzahl von Ausgleichsprozessen vermag man mit dem 
bisher von uns verwendeten Differentialgleichungssystem (3) zu erfassen. Dif- 
ferentialgleichung (4) entfallt ja, da wir — laut Definition — unter Ausgleichs- 
prozessen solche Umsetzungen verstehen wollen, die ohne augenblicklich ein- 
wirkende Strahlung vonstatten gehen. Wir nennen sie 


; 
j 
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a) Riickumwandlungen. Aus dem System (3) wird, da J = 0: 


es [aaa 
Ur 2 
\ 


: Misi N; 
; lar —%,4 + m,] NA | oe i (11) 


a =O, 
also ein System von linearen Differentialgleichungen, sogenannten «inhomo- 
-genen Diffusionsgleichungen», wie die einfache Transformation offenbart : 


“a 


—- —pyt 
Nee ee: 5 | 


i, Ae o 12 
(5; -— 4) ee ¥ fiir ‘aac | a 
| 0 | =O. 

___ b) Nebenrethen. Ein Spezialfall von Ausgleichsprozessen verdient, hier noch 
-erwahnt zu werden, namlich die spontanen Umwandlungen der betrachteten 
Teilchen — nach einer gewissen Inkubationszeit — in einen neuen, nicht bereits 
-durchlaufenen Zustand. 

_ Eine Erweiterung unserer Umwandlungsschemata (1, 2) illustriert dies am 
_besten: 
Hauptreihe Nebenreihen 
(x) + a > (B) = (B 0) > (B81) > (B82) ... 


| 


| (13) 
(o — 1) +4, > (@) = (@ 9). | 
Hauptreihe: (2) ++4¢—> (o): NAw); Me™. | (14) 


Nebenreihe: (7) > (tk): N,,(t, 4); Mins “in- 


Man beachte: ; N; erfasst den Anteil der Partikel des Zustandes (7), der in 
“den Zustand (i — 1) zuriickfallt, w;9 N; aber den Anteil der Partikel des Zu- 

standes (7), der in den ersten Zustand (7 1) der Nebenreihe tibergeht. Nebenrei- 
hen sind nur méglich, wenn 
| Mi< Bik: (15) 


Das Differentialgleichungssystem (11) wird also modifiziert : 
Misi Nisa 
(5 —2#,A + Mio) N= | 
| 0 | = (Dis 


Bs | | 
(57 — 4;,4 Hit) Nix = Bie-y Nig-y fir ksi. 


4 
a 1 
: ] 
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Nebenreihen kénnen natiirlich auch wahrend der Bestrahlung auftreten, die 
Gleichungen (3,4) sind dann entsprechend abzuandern; wir brauchen dies hier 
nicht zu explizieren. : 
Auf dem nun von uns dargelegten Wege lassen sich bereits eine Reihe fiir _ 
die Strahlenbiologie und Radiochemie wichtigen Verhaltensweisen von biolo- — 
gischen Objekten und Substanzen in guter Naherung beschreiben, beispiels- 
weise das Abklingen der Strahlenwirkungen mit der Entfernung vom Strahlen-_ 
herd, das heisst die Eindringtiefe und der Einflussbereich der Strahlung, dann : 
die kritische Bestrahlungsdauer, mit deren Uberschreiten tiefgreifende Schadi- 
gungen auftreten, dann Regenerationsméglichkeiten und deren Dauer, um nur _ 
einige fiir den Strahlenschutz wichtige Fragen zu nennen. : 


3. Integrationsverfahren 


Wahrend uns zur Integration von (11) die klassischen Verfahren der Theorie — 
der Diffusion bzw. Warmeleitung (wie Fourier-Methode, Laplace-Transforma- 
tion, Warmepole, Integralgleichungen u. a. [1]8) zur Verfiigung stehen, sind 
diese Verfahren fiir St6érungsprozesse ungeeignet, da nunmehr keine linearen 
Differentialgleichungssysteme auftreten, sondern nichtlineare zu ihrer Be- 
schreibung herangezogen werden miissen. 


a) Linearisierung. Wird aber neben der Voraussetzung, dass nur ein 
kleiner Bruchteil der den Ko6rper treffenden Strahlung absorbiert wird, der. 
Grossteil also den Kérper durchtritt, noch angenommen, dass nur kurzzeitig 
bestrahlt wird, das heisst, dass die mittlere Lebensdauer der Zustande lang ist 
im Vergleich zum Bestrahlungszeitraum, dann kann man das Differential-_ 
gleichungssystem (3,4) linearisieren. 

Wir transformieren dazu das System (3,4) mit 


Wee = renee Nae Mos | 17 
Nene” Sa eS Rales of 
Te, t) = Jo — Tw, t). (18) 


Nun kénnen wir, da M; < Ny; I < J, fiir die in (3,4) auftretenden Produkte 
N, J setzen: 


N,— M, = = , 
ree oe M)(Jy—1) » NJ -JyM, 
| M; Jy— 1) © Jo M; 


) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 260. 
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be 


Damit wird aus (3,4): 


0 
4 (5 - t, A+, Jo) Ma = Up Mp tre Nol, 


: (20) 
Gre 2, +; Jo+ ui) M; = Mi41 Migr +%- 1M; —1> t+ 0. 
a 
i) 

4f 4 y,NyJ = Jo fre Ma— SM), 1) 
e Be ae 4 


Der Zwischenzustand (i = £) ist damit auch erfasst, wenn man in (20) statt 
AV ;-1):-¢ = No — M, setzt. 

_ In einigen sehr einfachen Fallen gelingt eine geschlossene Integration des 
Systems (3,4). Allgemein jedoch muss man sich mit einer mathematischen 
Approximation der Lésungen begniigen. Die wachsende Verbreitung von 
-elektronischen Digitalrechnern legt natiirlich eine numerische Integration des 
gerade interessierenden Trefferproblems nahe. Der Nachteil einer solchen Vor- 
_gehensweise liegt oft darin, dass dabei der Uberblick iiber die funktionale Ab- 
hangigkeit der Lésung von den verschiedenen Parametern verlorengeht, der 
‘iibliche Einwand also. 

_ Eine Schwierigkeit ist — auch fiir den Digitalrechner — darin zu sehen, 
“dass ein Anfangs-Randwertproblem vorliegt. Der programmierungstechnische 
Aufwand und damit die Bearbeitungsdauer liesse sich aber in tragbaren 
-Grenzen halten, wenn es gelange, das Anfangs-Randwertproblem auf ein reines 
Anfangswertproblem zu reduzieren. Man miisste dazu zum Beispiel die gesuchte 
‘Lésung aus Funktionen aufbauen, die zwar das Differentialgleichungssystem 
‘und die Randbedingungen erfiillen, aber den Anfangswerten nicht zu geniigen 
brauchen. 

_ Im folgenden wollen wir nun — ausgehend von einer Bemerkung von 
-y.MisEs [2] — eine Approximationslésung fiir (3,4) konstruieren, die dann auch 
als Basis fiir die Programmierung fiir einen Digitalrechner verwendet werden 
kann, gewissermassen als «Unterprogramm», wie wir sehen werden. Die Ab- 
eitung geschieht der Ubersichtlichkeit wegen an einem einfachen Spezialfall ; 
‘die Verallgemeinerung der Vorgehensweise ist danach evident. 


__ b) Konstruktion einer Approximationslosung. Eine zwischen parallelen im- 
-permeablen Wanden eingeschlossene Lésungsmittelschicht (Figur 2) enthalte 
eine strahlungsempfindliche Substanz. Bevor auf das eingeschlossene Medium 
ein paralleler konstanter Quantenstrom tiberall gleicher Flachendichte J, trifft 
(Parallelstrahler), mag N,(x,t) = Ny > 0, die lokale Konzentration der strah- 
lungsempfindlichen Substanz im Pe aesuetand, raumlich und zeitlich kon- 
stant sein und die Konzentration N,(x, t) = 0 fiir ¢ < 0. Der Bestrahlungszeit- 
‘raum sei durch 0 < ¢ < ¢, gegeben. 


We 
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X=0 sail 


Figur 2 
Strahlung der Intensitat Jy trifft auf eine strahlungsempfindliche Parellelschicht der Dickel. — 


Mit der Transformation 
EootpeaN oie eae IN Gee 3 


Z(€,0) = MED; Fe,9 = Jed, 


Cake = a hash, 


reduzieren wir fiir diesen Fall im Differentialgleichungssystem (3,4) die Anzahl 
der freien Parameter: x,, Hg, Vos Yps Jy. No, taut x2, Te Th, 


0254. 0Zp _ ,, OZ Deen, 3 
ae aT Fea S477, Da Fiz. tv). @ 


Anfangsbedingungen: Z,(&, t) = 1; Z,(§,t)=0 fir t=0; 


Randbedingungen: JE s)= J, fir -P=0; 


OL. OZp __ ey 
athe) Oe tin = Se 0 aioe 


Fiir den Sonderfall x = 1 erhalten wir mit Z = Z, + Z p aus dem System: 


WA A 


We 0 Og 
mit Z = 1 fir t = 0 und 0z/dé = 0 fiir € = 0: L. 
Die Lésung dieser Differentialgleichung kann bei den angegebenen An- 
fangs- und Randbedingungen nur die triviale sein4) 


Zee en Se (23) 


4) Gotpicidhe [1], Seite 578. 
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Es ist also statt (22) nur das System: 


ea) a) 


OZ OPZ x oy a 
ea Se 2,J, Ge-—-J bt -9) 2] (24) 
m integrieren. Z, findet man dann aus (23). 

Ohne weiteres lasst sich die 2. Differentialgleichung aus (24) einmal formal 
‘nach € integrieren mit dem Ergebnis: 


a 


- 


" 
s 


T(E, 0) = Jpexp . pees) He Z,(é, 1) a (25) 


Die aus (24, 25) fliessende Integrodifferentialgleichung: 


; é 
0 x 5 a 7 " 
ee a =~ fy2Z009| v8 0-9) [216.9 a =—8",7) (26) 
; 0 


inearisieren wir dadurch, dass wir den auf der rechten Seite stehenden, an 
sich glatten g”(E, t)-Verlauf hinsichtlich + durch einen Polygonzug ersetzen 
(Figur 3). Bei t = 0 beginnend, errechnen wir 


gl (é) =e"(E,t% = 0) = Joe 


& 


Figur 3 

Treppenpolygon approximiert den glatten Funktionsverlauf. 
Wir verwenden diese — nur von é abhangige Funktion — approximativ fiir das 
_gesamte Zeitintervall 0 <1 <1, und erhalten so statt (25) die inhomogene 


 Diffusionsgleichung: 
j OZ a, 0°Z 5 ” 
a = —8(6) - 


ir 08 


; 
: 
: 


EDa es sich bei dem Integral dieser Differentialgleichungen um eine Approximation 
von Z, handelt, nicht um Z, selbst, fiihren wir ein neues Symbol ein: Z, > m®). 


(0) 02 (0) . 
an dia = en 
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Die Integration dieser Differentialgleichung geschieht mittels Fourier-Reihe, 
Zuvor transformieren wir (27) zweckmassig mit: 


MOE, x) = wO(E, 2) + gf) + ESF 2* goo) — FEE g(t) (28) 


auf die Form 


aw _ ow fe 0) 
dr oe? | 

mit der Anfangsbedingung: 
WOE, 0) =1—g(@) - EH) poy + FF gi(L) =O(€) (30). 


und den Randbedingungen: 


dw) a ‘ 
ye = 0) fiir é =e 0; ik . 


Die fiir das Zeitintervall 0 <7 <7, giiltige Approximation lautet nun in 
Fourier-Darstellung: 


wO(E, t) 


L 
L fase )ae +2 Ser Qainie cos F*. [ Oé) cos 18 gel (3 : 
0 


J 


Hieraus finden wir ein m(&, 7) = m(€). Dieselbe Vorgehensweise wenden 
wir an, um zu einer Naherungslésung fiir das Zeitintervall 7; <1 <1, zu 
gelangen, usw. 

Ersetzt man in (28) das (0) > (m) und den Index 0 >, so erhalt man 
zusammen mit: 


L 
wm(E, 1) = = hy Jae &) dé + 2 Se faiEt* cos PAS ['@,(8) cos LE a. 
eh . (32) 
@,(E) = mone) — g,(@) - ESE 42% pry y 2 te riz) , : 


giiltig fir t, St <t,,1, sukzessive die Naherungen bis zu einem gewiinschten 
Zeitpunkt ,. 


Diese Naherungslésungen bauen wir nun zu einem Kurvenzug zusammen: 


1) & Sel) — eltyaa)} mE, 2). (33) 


ay SE ae 
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-€(t) symbolisiert dabei die elementare Heavisidesche Sprungfunktion: 


7 rea oe 
E o(t) = ible Gee 

4 |o |<o 
‘ 


Eine gewisse Schwierigkeit liegt bei diesem Verfahren in der zweimaligen 
Quadratur der Funktion g/() > g,(€) . 

5 Man umgeht sie und ersetzt zugleich die Reihendarstellung durch eine 
Integraldarstellung der Approximationslésung, wenn anstelle des Differential- 
-quotienten 0m™/0r der zugehdrige linksseitige Differenzenquotient: 


| Mer, + _ mM, () (34) 
tritt. Aus (27) wird so 
1 a? n n 
<Ttt — Bm y1 = — 2 m, + nlf) = — h(E) » (35) 
R= 1/Ar. Die Einarbeitung der Randbedingungen bringt: 
; L 
my sa(6) = | G(E, 8) hyls) 4 (36) 
3 2 0 
_wobei G(é, s) die Greensche Funktion fiir die vorgegebenen Randbedingungen 
-ist: 
= Cosd (L — &) Cosas 

Vit) 5 OSss6, 

BS Kcash (Eas) Cos § ee 
OS — OS 
A She cSssl. 


Diese Lésung (36) ist auch fiir eine numerische Auswertung gut geeignet. 

Die Méglichkeit zur Verallgemeinerung zeigen wir am geometrisch gleich- 
_geformten Beispiel auf. Die Abweichung gegeniiber dem eben durchgerechneten 
liegt darin, dass nun % = %,/x, + 1 ist. Hier gilt (23) nicht mehr. Wir miissen 
also jetzt das System: 


Gea 8 o8z = : ~ y 
- 0&2 ae Zoo] = [(ea+¥ Zp a = —x#*§F (&, 7) , | 


: (37) 


) OZ, } 0?Z~ or fel, ra wo] [ea ¥ 2) a = +2 g"(6, t) | 
0 


Ot og 


“unter den oben genannten Anfangs- und Randbedingungen integrieren. g’ (&,7) 
_approximieren wir wieder — wie oben — durch ein Treppenpolygon und separie- 
“ren und linearisieren so das System (37) fiir das Intervall t, St S T+1- Jede 


' 
+ 
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| 
der beiden Integrodifferentialgleichungen von (37) ist damit unabhangig von. 
der anderen zuniachst fiir 0 <7 <7, lésbar. Am Intervallende 1, errechnen 
wir dann wieder die Funktionswerte m!?)(&, t,) ; m')(E, tT.) und nehmen diese als 
Anfangswerte fiir das Intervall t, < t S 75, usw. | 
_ Wir sehen, dass eine Verallgemeinerung auf mehrere Zwischenzustande am 
formalen Ablauf des dargelegten Weges nichts andert; die Integrodifferential- 
gleichungssysteme werden auf diese Weise stets fiir jegliches t, St <T,41 
linearisiert und separiert. 
Auch die Beriicksichtigung von riicklaufigen Spontanreaktionen, von Spal-_ 
tungen und anderem mehr bringt keine wesentliche Abweichung. (37) ist also” 
auch die Lésung fiir ein beliebiges Differentialgleichungssystem (3,4), wenn 
man h,(s) geeignet interpretiert. | 
Programmierungstechnisch gesehen, wird man also in allen durch (3,4) er-. 
fassbaren Fallen ein Unterprogramm fiir (37) benétigen; iiber h,,(s) wird dabei. 
zunachst natiirlich noch nicht verfiigt; G(&, s) richtet sich nach den Randbe- 
dingungen. Ein weiteres Unterprogramm ist fiir (33) herzustellen, auch das ist 
stets giiltig. Mit der jeweiligen physikalischen Problemstellung variiert bei 
gleichbleibender geometrischer Konfiguration allein das Programm fiir h,(s). 
Als weiterer Vorteil kommt noch hinzu, dass die Exponential- und Hyperbel- 
funktionen, wie sie in (37) auftreten, bei jedem Digitalrechner als «Archiv- 
programm» niedergelegt sind. 


@ Strahler 


O Medium 


a) 


Figur 4 


a) Parallelschicht, b) achsensymmetrische, c) kugelsymmetrische Einbettung des Strahlers in das 
strahlungsempfindliche Medium. 


Bei Koérpern beliebiger geometrischer Gestalt treten nach einem gegentiber 
dem obigen etwas abzuwandelnden Lésungsverfahren statt (37) Mehrfach- 
integrale auf. Man wird also versuchen, soweit es die Umstinde irgendwie 
gestatten, mit den geometrischen Grundkonfigwrationen (Figur 4) auszukom- 
men, um so nur eime Ortsvariable in der Rechnung beriicksichtigen zu miissen. 

Im Falle, dass ein achsen- oder punktsymmetrischer Strahler ausserhalb 
des zu betrachtenden Mediums liegt, kann man fiir einen nicht allzu grossen 
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: 

Winkelbereich gy in erster Naherung folgende Anordnung durch die Grund- 
_konfiguration (a) ersetzen. 


a el 


Figur 5 
Blende von kleinem Offnungswinkel m erzeugt approximativ eine Parallelstrahlung. 


Bevor wir nun zu Beispielen iibergehen, sind noch die beiden folgenden 
Fragen zu beantworten: 
; 1. Existieren zum System (3,4) in0 <*« </;0 <t <#, stetige und stetig 
- differenzierbare Lésungen mit den gewiinschten Anfangs- und Randbedin- 
_gungen ! ? 
; 2. Konvergieren unsere Naherungslésungen gleichmassig in « und y mit 
_ At > 0 gegen die Lésung von (3,4) ? 
| Man kann zeigen, dass die Existenz von Lésungen der gesuchten Art ge- 
sichert ist [4, 5]. Wir verzichten hier auf die Wiedergabe eines Beweises. 

Die Frage 2 beantworten wir mit dem Beweis des 

Satzes: Der Unterschied U, +4(%) = N(«,t,41) — 0(%, 441) = Ny4i.— M41 

geht fiir» = 0; 1; ...; [¢,/Af] mit At gleichmassig in x und v gegen Null. 

Bewets: Der Bestrahlungszeitraum seid <<t<t.. N, g(x, N), 0g) /ON seien 
_stetige und mindestens einmal stetig differenzierbare Funktionen ives Argu- 
‘mente in0<*<1;0<t<4. U(x) =0. 

Setzen wir in die Differentialgleichung 


| oN ON ON) | Noa — Ny _ 
oer [a ag oe 8) 


, ein, dann geht bei stetigem N der Lim 6, 41(%) > 0, und es ist zum Zeitpunkt 
¢=t,,, mit A? = 1/A¢: od 


@Ny41 _ 92 (N 


dx? y+1l i, _ go 41%) => 6, 41(%) = 


Wir subtrahieren davon nach (35): 


d? Ny. N 
ayes - #? (My 44 y. n,) af g, (x) 
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BU y41 _ ge U, 41 = —A U, + [e"(%, Ny 41) — 8"(%, %)] + O,41(%) 


di 39) 
= —H,(x,). | 


und erhalten: : | 
| 


U, +1(*) geniigt also einer Differentialgleichung vom Typ (35). Die zugehérigen 
Randbedingungen sind dU, ,,/dx = 0 fiir x = 0; 1. : 
Um die gemass (36) gefundene Lésung 


1 
Upsalt) = 3 | Gle, 8) HE, A) ae (40) 
0 
abschatzen zu kénnen, fiihren wir die Hilfsfunktion V(x) ein, iiber 
V'—-#BV=—|H,(x,)| und 2° —0 far x=0; 1. (41) 
Nach (36) ist dann 
1 
ae, wie 
Vix) = Ge f Gle, 8) |HE a)| df SU, ale). (42), 


0 


Man erkennt aus (41), dass V(x) > 0in0 <x <J, daanallen Wendepunkten 
der Hilfsfunktion V(x) > 0 sein muss und an dazwischenliegenden Minima 
(das heisst V" > 0) fiir V(x) nur positive Werte eintreten konnen. Die Rand- 
werte verlangen an den Intervallenden horizontale Tangenten. Es liegen dort 
also entweder Extrema oder Wendepunkte vor. Also ist in 0 < x </ stets 
V(x) 20. An der Intervallstelle, wo | V(x) | einen Maximalwert annimmt, ist — 
| (x, A) | 


Fig) | eae 


| H,(x, a) | — 22 V(x) = —V"(x) > 0 > 


Die damit gewonnene Abschatzung fiir U,41(%) lautet dann: 


| U,4a(0) | S sy Max|H,(x, a). (43) 


Lasst sich noch zeigen, dass H y(x, A) fiir wachsendes 4 beschrinkt bleibt, dann © 
ist der Beweis fertig. 


Wir setzen in (39) nach dem Mittelwertsatz: 


W W 0 v 
a", Ny +1) <8", »,) = (N41 — ,) So, (44) 


0g" /ON ist die Ableitung an einer zwischen N y+1 und n, liegenden Stelle, und 
erhalten: 


H,(x, 4) = # U, — (N,44— N, 4+ U) oe — 6, 41(%) . 


A 
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~ 


Da die auftretenden Funktionen im betrachteten Intervall stetig und minde- 
-stens einmal stetig differenzierbar sind, das heisst 

"| dg" é A[/ oN A 

Bfan| <4: [Poe | Ae oe (GP) = 15, 41(%)| S52 [Bt %0,41(%)] 


by 4d 


“mit konstantem A; B, unabhangig von »; 4, gilt mit U,(x) =0; »< Az: 


t 


M+) | S (1 i +) Max | U,(x) | (mod =) 
< (1+ 4) Max|U,(x)| (mod {). 


c) Lésung bei Ausgleichsprozessen. Aus den oben aufgezahlten Methoden, 
‘die zur Integration des linearen Differenzialgleichungssystems (11) zur Verfii- 
-gung stehen, wahlen wir den Fourier-Ansatz, weil er — im Hinblick auf die 
_praktische Auswertung konkreter Beispiele — der einfachste ist. 

_ Zweckmiassigerweise beginnt man mit der Integration der Differentialglei- 
chung fiir den Endzustand (w). Sie ist ja dem System nicht angekoppelt. Die 
' daraus errechnete Lésung setzen wir in die Differentialgleichung fiir den Zwi- 
schenzustand (w — 1) ein, usw. 

; Fiir die Grundkonfiguration a (Figur 4) wollen wir die Vorgehensweise dar- 
legen. Zu integrieren ist hier zuerst die Differentialgleichung 


ON O'My 0 ; ( 45) 


nn oe 
“mit On,/dx =0 fir x=0,1 und n,(x, 1) =©,(x) fair ¢ = 0. 


Wir setzen — gemiss den Randbedingungen — an: 


nip >? ed) aay ae 
k=0 


und finden so aus der Differentialgleichung (45) und der Anfangsbedingung: 


} 1 
. 1 
AW) = 7 [O) a 
i 0 
1 

2 —(2km/l)*x x 

AVP(t) = 3 e~PAaiD ot [O@ cos ARS de, kh +0. 

z; 0 
-Damit gehen wir in die Differentialgleichung: 
2ku%x 


ONw_ OP Ny_ — t w) 
5 a HX y-1 ae = Fo eWHo—t Ha) DA (t) OED aaa ; 


| Die Integration dieser Differentialgleichung geschieht nun wiederum mit einem 
die Randbedingungen identisch befriedigenden Fourier-Ansatz, usw. 
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Fiir den z-ten Zustand gilt die Differentialgleichung 


i. 2. eo Ae a ; 2k ex 4 


Die Lésung dieser Differentialgleichung ist bei der Anfangsbedingung 


n,(2e, 0) = 3,'al? cos *42* (47) 


k=0 


gegeben durch: 


1;(X, t) = DAI (2) as ah 
ko 
. (48) 
a + Hiss [elitr De eae hats ACLS are 
6 


Ali) (t) =e @halb xt | 
Da alle A((t) einfache Exponentialfunktionen sind, bietet eine solche rekur-- 
sive Integration des Differentialgleichungssystems (11) keinerlei Schwierig- 
keiten. Auf die Ausfiihrung anderer Integrationsmethoden wollen wir hier. 
verzichten. 

Bei achsensymmetrischen bzw. punktsymmetrischen Anordnungen treten 
anstelle der Winkelfunktionen Zylinder- bzw. Kugelfunktionen. Die von uns 
fiir den Parallelstrahler abgeleiteten Integrationsverfahren sind nach einer 
geringfiigigen Abanderung auch fiir diese Falle gut brauchbar. 


d) Beispiel. Vorgelegt sei die Anordnung von Abschnitt 3b) (Figur 2). Die | 
Diffusionsgeschwindigkeiten dieses Eintreffervorganges seien einander gleich, 
x = 1, die Teilchen des Zustandes (8) mégen eine gewisse Inkubationszeit 1, 
bendtigen, bevor sie sich in den stabilen Endzustand (8 1) verwandeln. Der 
Bestrahlungszeitraum falle mit der Inkubationszeit zusammen. 

Die zugehérigen Gleichungen sind fiir: 0 < + < ..: 


aZe 07Z ; : Ps 
Z,+ 2Z,=1, eT. ~hZp0| 8-00) [2,018 (49) 
0 


fir: Z,= 1, = 0; 0Z,/0§ =0 fir €=0; L. 
Bir. 4° nd Hpi = Hy X% p13 My = X,(¥_ No)? [Lp: 


io; Ope ee 02h, ee. On eke 
Ooo peiuas ee Zp» ie “paper = Be Zp (50) 


mitZ,=@(€); Z,,=0 fir +=7,; 0Z,/0§ = 0Z,,/0§ =0 fir €=0; L. 
Mittels numerischer Integration von (49) nach obigem Verfahren finden wir LF 
Dieser zahlenmiassig vorgegebene Funktionsverlauf wird der tiblichen Fourier- 


~ 
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Analyse unterworfen, mit dem Ergebnis: 


; P(E) = 3) a, cos 4 ae 5 (51) 
: R=0 


Er dient uns danach als Anfangswert fiir das System (50). Man verifiziert 


durch Einsetzen in (50), dass firt, St St, 


Z g(t) = D>, eres ar ee) cos 


k=0 


2kaé 
et 


(52) 


-Damit errechnen wir aus der 2. Differentialgleichung (50) unter Berticksich- 


tigung des Anfangswertes 


Zoslé, 2) — 3/ A,(t) cos 225 , | 
k=0 
— Hp ak ‘ - Fp 1 (4h n*/L*) (t—-7,) __ p—- (4 Rre/L*) (c—ty 
Ax(e) Hea — (%py — 1) 4k? n?/L? pe pence tap bre 


‘Diese Fourier-Koeffizienten konvergieren, wie man leicht iiberpriift, minde- 


stens wie k~4, also sehr gut. 


Abschliessende Bemerkungen und Zusammenfassung 


Mit Hilfe der «Treffertheorie» gelangt man bei einer grossen Anzahl prak- 


_tisch wichtiger radiochemischer Umwandlungsvorgange zu quantitativen Aus- 


sagen. Wir zogen deshalb diese Theorie zur Bearbeitung von strahlenbiologi- 
schen Prozessen, wie sie neuerdings ganz besonders in der Raumfahrtmedizin 
auftreten, heran. Natiirlich sind die durchgefiihrten Betrachtungen keinesfalls 


allein auf die Astronautik beschrankt. Schon bei einfachen rontgenologischen 


“= 


Untersuchungen oder solchen von radioaktiven Einwirkungen ist die Proble- 
matik ahnlich. 

Das mathematische Erfassen radiochemischer Umwandlungen kann mit 
Hilfe partieller Differentialgleichungssysteme geschehen, deren Herleitung — 


_ peschrankt man sich auf das Prinzipielle und abstrahiert von Nebeneffekten — 


_ keine allzu grosse Miihe bereitet, deren Integration aber, nahezu in allen Fallen, 


nicht mehr geschlossen zu bewaltigen ist, also approximativ durchgefihrt 


werden muss. Wahrend bei den sogenannten Sekundarreaktionen (Ausgleichs- 


prozesse) die beschreibenden Differentialgleichungssysteme linear bleiben und 
so mit den iiblichen mathematischen Lésungsverfahren bearbeitet werden 


_ kénnen, verlangt die Integration bei den sogenannten Primarumwandlungen 


_ (Stérprozesse) einen ausserordentlich umfangreichen mathematischen Aufwand, 


da hier die Differentialgleichungssysteme nicht mehr linear sind. Man ist in der 


- Regel, will man irgendwelche Aussagen fiir den Reaktionsablauf erhalten, auf 


(53) 
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eine numerische Integration angewiesen. Die gebrauchlichen Differenzenver- 
fahren sind hier jedoch ungeeignet, wie man leicht zeigt; jeder Kompromiss 
zwischen Genauigkeit der numerischen Ergebnisse und Bearbeitungszeit, selbst 
beim Einsatz elektronischer Digitalrechner, bleibt unbefriedigend. Eine nume-. 
rische Bearbeitung von Mehrtreffervorgangen bendtigt auf diese Weise ausser- 
ordentlich lange Rechenzeiten, handelt es sich dabei doch um Randwertauf- 
gaben, bei denen mehrere Parameter durchvariiert werden miissen, um die 
Randwerte zu erfiillen. } 

Um einen Uberblick iiber den Parametereinfluss zu erhalten und um diese : 
Jangen Rechenzeiten abzukiirzen, wurde ein geeignetes Integrationsverfahren 
entwickelt. Die Uberlegenheit dieses Verfahrens hinsichtlich Genauigkeit und 
Schnelligkeit zeigt sich ganz besonders bei den Mehrtreffervorgangen. Die - 
Vorgehensweise dabei ist der Ubersichtlichkeit wegen an geometrisch einfach 
geformten Anordnungen hergeleitet und illustriert. Einer Verallgemeinerung — 
auf kompliziertere Probleme steht prinzipiell nichts im Wege, da das Integra- | 
tionsverfahren sehr flexibel ist und den Rand- und Anfangsbedingungen einer 
grossen Anzahl von Situationen angepasst werden kann. 

In einer folgenden Untersuchung werden damit konkrete Aufgaben bear- 
beitet und eine Fehlerabschatzung fiir die Ergebnisse angegeben. 
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Summary 


Many substances, when subjected to radioactive radiations, undergo definite 
changes. With the help of the Target-Theory it is possible to discuss some special — 
effects such as the depth of penetration of the rays, the allowable doses of radiation 
energy, the after effects, protection against radiation and so on. This paper deals 
with interference processes, i. e. processes where the radiations are incident; and 
the compensation process, i. e. processes that follow afterwards. The solution of the 
corresponding differential equations by means of a digital computer is outlined, 


(Eingegangen: 29. September 1960.) 
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“Magneto-Fluid-Dynamics of Thin Bodies in Oblique Fields 


, By KeitH Stewartson, Durham, Great Britain’) 
a, 


1. Introduction 


_ This paper is concerned with the steady motion of a fluid which is incom- 
pressible, inviscid and electrically a perfect conductor, past a thin non-conduct- 
ing body in the presence of a magnetic field which when undisturbed is uniform 
and orientated in an arbitrary direction. The interest in this problem is partly 

intrinsic and partly because it illustrates the effect of obstacles and boundaries 
_on the general motion of the fluid. In addition it raises questions of the correct 

posing of boundary conditions and of the interplay between main stream and 
boundary layer which at the present time cannot be settled by reference to 
experiment but only by ensuring the complete consistency of the assumptions 
made. 

The velocity and magnetic fields have been determined by SEARS and 

_ RESLER [3]2) if the undisturbed fields are either parallel or perpendicular 
‘while an alternative view of the solution in the case of parallel undisturbed 
fields has been put forward by STEWARTSON [4]. The object of the present paper 
is to throw further light on the flow pattern, by considering a more general 
_ problem of which the problems above are special cases. 
In general the boundary conditions specifying the problem contain an 
arbitrary constant (STEWARTSON [4]). This constant is however known for a 
real fluid of high electrical conductivity o and low kinematic viscosity » and 

if our fluid is regarded as the limit of a real fluid then the arbitrary constant 

is equal to 

Lim 42 (vo o)t# =A 


o— 00 
v—>0 


where o is the density of the fluid. Now 1 is the ratio between the discontinuity 
in the tangential components of the magnetic field at the surface of the body 
to the tangential velocity of the fluid there. Further, in the theory it is assumed 
that all disturbances caused by the body are relatively small and hence it 

follows that unless 4 is vanishingly small this assumption breaks down. We 
conclude that for a general prefectly conducting inviscid fluid the disturbances 
_ in the steady state motion of the fluid are not small, not only when the undis- 


1) Department of Mathematics, Durham Colleges in the University of Durham. 
2) Numbers in brackets refer to References, page 271. 


a 


262 Kr1rH STEWARTSON ai | 

| 
turbed fields are aligned, as has been suggested earlier, but also for oblique 
fields. There remains only the case 2 = 0 and fortunately this is the one which 
is most likely to be of interest from an experimental point of view. We shall 
show in this paper that if A = 0 the disturbances produced by the body are 
small in general, consisting of a potential field and two waves extending te 
infinity on either side of the body directions on the downstream side of that 
of the magnetic field. In particular the solution given by SEARS and RESLER [3] 
for a perpendicular field is recovered as a special case. However, the solution 
obtained does contain an arbitrary constant of a similar kind to that occurring 
in conventional aerofoild theory where it is removed by the Kutta condition, 
of finite velocity, at the trailing edge. This condition in aerofoil theory is in. 
agreement with experiment and supported by a considerable theoretical 
argument. Here there are no experiments available and a corresponding 
theoretical argument involves difficult questions of magnetohydrodynamic | 
boundary layers which are not yet resolved and are beyond the scope of this 
paper. We do show, however, that the application of a Kutta condition at the 
trailing edge leads to no difficulties if # > 1 where m is the ratio of the speeds 
of the fluid and of Alfvén waves. If m < 1 however, there is an infinite discon- 
tinuity in the forces on the body when the angle between the directions of the 
magnetic field and the velocity field is 1/2 arccosm? unless there is a local 
transference of the Kutta condition from the trailing to the leading edge of the - 
body. 

Apart from this case the solution breaks down only when the undisturbed 
magnetic and velocity fields are parallel, both wavelike components to the flow 
being formally infinite. As the angle between the magnetic and velocity fields, 
arctan f, tends to zero it is found that increasingly large disturbances are 
propagated to infinity without any diminution, approximately upstream and 
downstream of the body if m < 1 and downstream only if m > 1. Further, the 
drag on the body also tends to infinity. The theory is not completely valid 
when f is of the same order or less than the slope of the body but it is considered 
that it does provide support for the view of the solution, when the field is 
aligned, given by STEWARTSON [4]. 


2. Governing Equations 


Consider a perfectly conducting inviscid incompressible fluid of density 0 
in steady motion. Relative to fixed axes O x yz let the undisturbed velocity 
of the fluid be (U, 0,0) and the undisturbed magnetic field be (« H, 8 H, 0) 
where U, H, «, 8 are constants, B > 0 and a2 + B? =1,.A thin non-conducting 
cylindrical body is fixed in the fluid near O so that the equations of its upper 
and lower surfaces are 


y=f@), |Z]|Sa, y=f@), |x| <a, (2.1) 


NEN NE 
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Tespectively f(x) being small. The presence of the body disturbs both the 
magnetic and velocity fields; let us write 

a v= U(1+4+ 4, Uy, 0) (2.2) 
for the velocity field and 

4 H=H(a +h,, p +h,, 9) (2.3) 


for the magnetic field at a general point of space where u,, u,, 1, , are assumed 
‘to be small quantities whose squares may be neglected. With the same assump- 
tions as in STEWARTSON [4] MaxweELL’s equations become 

curl (v , H)=0, divH=0 (2.4) 


. 
4 


and EvLER’s equations, modified by the addition of the electromagnetic body 
force, 

4 gu 1 ; 

Ox curly = —grad(P eas v*) ~ Cxare curl H , H, divy=90 (2.5) 
where p is the pressure and the permeability is unity. Substituting from (2.2), 
(2.3) and neglecting squares and products of small quantities we get 


3 hy Oe OU, Ohy _ _ Oty Oy 
4 a a era ra ae ae | 

1 pee I eee (tae Oe) 

: Ox o ox m= \ dx dy }? 

4 (2.6) 
Ply oee 3 deny (2 s aha 

4 Ox o Oy m2 \ 0x dy }’ 

Ng Ohy Og Oy _ 

aor oy ve rae oy 

~ where 

f 2 

rer Pod (2.7) 


so that m is the ratio of the undisturbed fluid speed to the speed of Alfvén 
waves. ; 

In view of the symmetry between x and y in these equations we can expect 
that in a typical solution all dependent variables are proportional to F(% — ¢ y) 
where F is an arbitrary function and c is a constant to be found. Thus writing 


: h, = F(x — cy) 
_we find successively that 
+f 
| ae. owin. 
h=cr, ate rae By pgs 
| Pe fe we o:2 2)| F 
o U? Blac ont (1+ 0] 


nial one te eee be Shee | 
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whence from the last equation of (2.6) | 


7 


e+1l=0 or (Be-—ajt=m. (2.8) 


If c? + 1 = 0 all the dependent variables are harmonic while if B c — # = + m 
they are wave-like. Thus in the steady solution there could be two waves in 
the region y > 0 and two waves in the region y < 0. In order to decide which 
are relevant to the present investigation we must specify the problem more 


undisturbed 
velocity 


Figure 


Zones of wave-like solutions (shown shaded). 


precisely. Let us suppose that initially the body is moving with the fluid, 
with velocity (U, 0,0), and that the magnetic field is uniform. At time ¢ = 0 
the body is suddenly brought to rest, in the configuration described by (2.1). 
This problem is similar to that discussed in STEWARTSON [4] and a parallel 
mathematical argument can be used to determine the ultimate steady flow. 
There is no need to give it here: the conclusions are identical with those of the - 
following simple argument in physical terms. At the instant ¢ = 0, Alfvén waves 
are set up at the surface of the body and travel outwards, one in the direction | 
of y increasing and one in the direction of y decreasing. Were it not for the 
motion of the fluid they would travel along the magnetic lines of force but 
in fact they are deflected downstream by convection. Ultimately therefore the 
wave-like part of the solution will consist of two Waves one on each side of 
the body and on the downstream side of the lines of magnetic force crossing it. 
The corresponding values of ¢ are 


Oy ag fay 04 0. a ee oe Oe (2.9) 


. 
4 


is 
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The complete solution can be written as 


0d c 0 1 
= Gh ae cy), uy = 5 — J i GY), 
‘ Od F) 0 
bea St By tee Fe—e49), hy=—B So +a 5h + Ble cxy) 54 (2.10) 
i 
PP? 0 1+ 
0 i = - as ae F(% — c,/¥) 


if y > 0 where 


“and fo is the undisturbed pressure. If y < 0 yi c; by c_, Fby Po ane 
m by — m. 
3. Boundary Conditions 


On the body the normal component of velocity is zero and hence 


w=, (4) o y=0,, |x| Se (3.1) 


as is usual in slender body theory. Inside the body since it is an insulator 
divH=cnlH= 0. 


Accordingly if h, is to be small in the body 0h,/dy must also be small and so 
the change in h, from one side of the body to the other must be second order. 
Similarly the change in h, from one side of the body to the other is second 

order. At the surface of the body the normal component of the magnetic field 

“must be continuous and, since the normal component is non-zero, the discon- 

_tinuity in the magnetic field is proportional to the tangential component of the 
fluid velocity. In fact since the undisturbed magnetic field is everywhere 
constant we have 

Aah, mh, at ye One leh 8, | 


AU (l+u,) =Hh, , —Hh aty=0_, |x| Sa | 


“body *fluid 


(3.2) 


‘where A is discussed and defined in the introduction and the usual approxi- 
“mations of thin body theory have been made. The presence of the undisturbed 
fluid velocity in (3.2) means that in general the assumption of small disturb- 
ances fails. Fortunately in the flows which are likely to be of greatest interest 
‘Ais extremely small and so here we shall consider the special case A = 0 when 
‘boundary conditions require the magnetic field to be continuous at the body 
‘surface. Combining this with the earlier conclusion that the magnetic field 


a | 


| 
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inside the body is independent of y to first order it follows that the nidgnes 
field in the fluid must satisfy 


Lim h,=Limh,, Limh,=Limh,, if |x| Sa. (3.3) 


y>0+ y—0- y— 0+ y>0- : 

On the plane y = 0 but away from the body it follows from the fundamental 

equations that u,, u,, h,, h,, p are all continuous and so we conclude that : 
F, (x) = F_(x%) =0 


if | «| > a and that the singularities of ¢ are either inside the body (i. e. on” 
|x | <a, y=0 to the order of the approximation) or at infinity. Finally at 
infinity grad ¢ > O because all disturbances originate at the body. 


4. The General Solution 
Let 


P= b+ $, (4.19 


where gy is an odd function of y and ¢, is an even function of y and let 


0 0 
(fo = Agia) sgny, 0 — Bix), 8 = 4 (x), 8 — Bix) sgny, (4.2) 


when y = 0, | «| Sa. Then from the properties of harmonic functions 


= — | A(x’) tan-1 7—* dx’ 
o(% ¥) = [aut ie cer aa A 
(4.3) 
b(t, y) = Fe | Be ) log [( — x')® + 9) dx’, 
and in particular 
5 ode A Mele) dee ee Brae 
Bilas f SS, Aya | BNO SS ee 


In terms of Ay, A,, By, B,, F,, F_ the boundary conditions at the body 
become 
1 ' 
Bo(x) + B,(x) — — Fy (x) = f', (x) , 


m + 


— 


from (3.1) and (4.5) 


from (3.3). 


et 
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E A second relation, in addition to (4.4), between A, and By can now be 
found from (4.5) by elimination. It is 


Um (m? + 1 — 2.2) Bylx) — B (1 +m?) Adle) 

= m (m2 + 1 — 202) P(x) — a m? Q(x) = R(x) a 
D here 

| 2 P(x) =f, (0) +f), 200) =f.) -— 7). eee 


The solution of (4.6) subject to (4.4) is known apart from an arbitrary 
constant (Rorr and CHENG [2)) and is 


a 


rR real | ae) 


z Ls a+ a— x’ 

“q —a (4. 8) 
, ae Waa aCsin@ x 

a te ae ~ wi-d 0 

A (a — x) (a + *) 

where Cis an unknown constant and 

® Pg eee he oe Ga: (4.9) 


B+ mm)’ 


: It is natural to fix C by means of a Kutta condition at the trailing edge 
x = a, implying that C = 0, but this means that there is a violent change in 
“the solution as m? + 1 — 2? changes sign since 6 jumps from 0 to 1. 
Another possibility is to assume that the Kutta condition is applied at the 
trailing edge if m+ 1 > 2a? and at the leading edge if m2 +1 < 2 a? which 
“is equivalent to requiring @ to lie in the range (— 1/2, 1/2) and C = 0. The 
“precise condition on C can only be found by considering the limit o > og, 
_y > 0 and in particular the boundary layer at the body. Such an investigation 
is difficult as indicated by the work of CARRIER and GREENSPAN [1] and must 
“be deferred for the present. It is of interest however to note that if o v= co 
it is easy to infer from their paper that if a flat plate is moving lengthwise 
“with velocity U parallel to a magnetic field of undisturbed strength H then 
(i) if m > 1 the boundary layer on the plate starts at the leading edge and 
ends at the leading edge, 
d (ii) if # < 1 the boundary layer on the plate starts at the trailing edge and 
ends at the leading edge. 
~ These conclusions support the view that C = 0 and WO 2, “but it must 
be remembered that in the flows described in this paper ov = 0. Further as 
we shall see below the argument of the present paper breaks down if the 
undisturbed fields are parallel. 


, 


# 
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Once A,(x) has been found the other unknowns follow from 


(m2 +1 — 2a?) F(x) =m (20%? -—am—1) Q(x) + BI—2am-+ m?*) A(x) 
(4.10) 


(m2 + 1 — 2a?) F(x) =m (202 +am—1) Q(x) —B(14+2am-4 m?) Aol 


2 ae | 
| 
(m? + 1 — 2 a2) B(x) = m2 Q(x) — 2a B Ag(x) | 

The Ade field is now known (apart from C) and the lift and drag form on 
any body may now be worked out. Thus on a body of zero thickness, for which | 
O@) = 0, the lift 4 


(mm? + 1)* — 4 a? mt 


L=— [ (by-04 — Pyao-) dx = 20 0 & me gna Ra [Ae ) dx (4.13), 


—a 


and the drag 


D= | (by-04 fi, — by-o- tL) dx | 
ry (4.14) 
ey 6 ucts Nh [a0 


OG oo 


In particular for a flat plate at incidence ¢ 


Pifee SAT ee and a—x\é aCsin@ a (4.15 
(x) é o(x) = esindz (- r =) Ne TP (4.15) 
T Ms (m? + 1)2 — 4.02 m? 
Few 8 sak ee ge Pee) (ty 


and D=eL. Since 6-0 from above as m?+ 1-22 from above it fol- 
lows that to keep L, D finite C >0 as m2 + 1-—->2.«2 from above. However 
9 ->1 from below as m? + 1-> 2? from below and so to keep L, D finite 
C>—2Ze as m*+1->2«* from below; from (4.15) this suggests that the 
Kutta condition is transferred from the ee edge to the leading edge as 
m* + 1 decreases through 2, or at least that the stronger singularity in 
velocity is transferred from x = a to x = — a. 

Again for a symmetrical body in which P(x) =0 and f,(x ) = file 
x = + a, the lift is still given by (4.13) while the drag 


nee Us tx a mm [ (mm? + 1)* — 4 a8 m3] 
iO Bm (m? + 1) [oe is mi+1— 2 a [pee eo | 
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the corresponding expression in terms of A, being somewhat more complicated. 
In particular for a body consisting of two circular arcs and of maximum 
thickness ¢ a, Q(x) = — ex/a, 


___ € *+2a6)amsinn@ (a—x\8 ail ean 
a, a? m+1—2a2 (-—*) (pea? (aaa? (4.17) 
/ ’ 
1, (m? + 1)? — 402 m? O Ot eet 
20Wxa m? (m? + 1 — 2 a) | erp (4.18) 
Bm (m* +1) D a 2 a(m? + 1)2 — 4 a? m?) 
20 Ua 3 z (m2 +1 — 222)? 
ak 2 92 Le gel = ‘4 
x [= om? (1 ae ON’ ) + Bex (1 +m) C (1 2 6)]. 


In order to preserve finite values for L, D in this case as m?+1>2«? 
from above it is clear that C (m? +1 — 2.«?)-* must remain finite. On the 
other hand as m?+ 1—>2.«? from below, 


Zeam 
e m+ 1— 2a? Au 

since 9 > 1 from below. Again this suggests that as m? + 1 decreases through 
the value 2a? the stronger singularity in velocity is transferred from the 
trailing edge to the leading edge. 

It is noted that if the Kutta condition is transferred to the leading edge 
‘when m2+1< 2a? then the drag given by (4.19) is always positive since 
}6| <1/2andC=0. 


5. Magnetic Field almost Aligned with the Stream 


This special case is of interest as it throws some light on the motion engend- 
ered by a body moving in the same direction as the imposed magnetic field, 
about which conflicting theories have been put forward. An aligned field is 
‘obtained by setting 6 = 0; taking f small, « ~ 1, the present theory gives 


m (m2 — 1) By— B (1 + m?) Ay = m (m? — 1) P(x) — m2 Q(x), (5.1) 


so that either the potential contribution to the solution remains O(e), where 
-P, Q are both O(e), or its dominating part is the homogeneous solution of (5.1) 
A, = 4 C/(a? — x*)4? with C = O(1). While this is a mathematical possibility 
it makes no material difference to our qualitative conclusions below and so we 
shall exclude it and suppose Ay, By are both O(c). Accordingly when f is small 


F (x) = _m O(#) F,(x) = — m Qi) (5. 2) 


= i a 1+m 
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other contributions being O(8). Thus a body without thickness gives rise to 
velocity and magnetic fields which are formally O(e) in the limit KC 0. How- 
ever, any body with thickness, so that Q + 0, gives rise to a velocity field in 
y= 0 OL 


n= gmt OG), 4 = Ol) 6.3), 
and a magnetic field 
b= — Fry Oe), hy = Ol) (5.4) 
Further the pressure on both sides of the body is 
by + “2 Ola (5.5) 
leading to a drag . . 
De 2m [ Q2(2) dx (5.6) 


which + oo as B > 0. 

The solution suggests therefore that as the fluid becomes aligned the 
perturbations remain small except within two bands of widths 2 a B/(1 + m) 
originating at the body. If m < 1 one extends upstream and one downstream 
of the body while if m > 1 they both extend downstream. In these bands the 
pressure and the components, parallel to the undisturbed direction of flow, 
of both the magnetic and velocity fields are large, the other components 
remaining small. 

These conclusions are in agreement with a theory of the motion of a body | 
in an aligned field put forward by the author STEWARTSON [4] in which it was 
argued that the velocity and magnetic fields cannot be only slightly disturbed 
everywhere and that directly downstream of the body there would by a finite 
disturbance to the fields even at a distance from it; there would also be a 
finite disturbance at an infinite distance upstream as well if m < 1. 

They appear to be at variance with an alternative theory put forward 
by Sears and RESLER [3] in which they argued that the velocity field would 
be only slightly disturbed everywhere and identical with that when H — 0, 
while the magnetic field would be parallel to the velocity field outside the body 
and zero inside the body. It must be remembered however that the theory 
when f is small is only valid if ¢/B is also small and cannot be wholly correct 
if ¢/B is large because the waves from the body would recross it at some stage. 
It is possible that h,, w, are complicated functions of e/B which are O(e/8) when 
é/B is small but + 0 as ¢/8 + oo. Intuitively such behaviour seems unlikely, 
raising further difficulties of a physical nature and it is considered that this 
solution provides further evidence in support of the author’s own theory. 


oy Aa\w 
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=, ; Zusammenfassung 


Es wird die stationare Bewegung eines vollkommen leitenden reibungslosen 
_Mediums um einen schlanken K6rper in einem schraggerichteten Magnetfeld be- 
-handelt. Es wird hervorgehoben, dass im allgemeinen eine der Grenzbedingungen 
zur Folge hat, dass die Stérungen durch den KGrper nicht klein sind. — In einem 
' speziellen Falle, der von grossem Interesse sein diirfte, gibt es keine «prima facie» 
_Einwande gegen die Annahme vom Vorhandensein kleiner Stérungen, und Losun- 
gen kénnen gefunden werden. Diese enthalten eine allgemeine Konstante ahnlich 
_derjenigen, die in der konventionellen Theorie der diinnen K6rper (ohne Magnet- 

feld) vorkommt und die in analoger Weise eliminiert werden kann. — Zwei spezielle 
_ Faille werden diskutiert, von denen der eine von Bedeutung sein kénnte fiir die 
- Theorie der Bewegung in parallelen Feldern. 


_ (Received: October 3, 1960.) 
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11. Sitzung des Schweizer Komitees fiir Optik (Sektion Lichtoptik ) 
vom 26. Oktober 1960 in Ziirich 


Am 26. Oktober 1960 fand im Physikgebaude der ETH in Ziirich die 11. Sitzung 
des Schweizer Komitees fiir Optik (Sektion Lichtoptik) statt. Anlasslich dieser 
_ Sitzung orientierte der Prasident iiber die Tagung des Internationalen Optischen 
_ Komitees vom September 1959 in Stockholm. 


Es folgten 3 wissenschaftliche Mitteilungen, von welchen nachfolgend eine 
_Zuasammenstellung gegeben wird. 


: MOd¢glichkeiten und Grenzen von Sehleistungspriifungen in Betrieben 
_ (mit Demonstration eines Versuchsinstrumentes der SBB). Von H. RUetscut, Basel’). 

Messungen von Sehleistungen an grésseren Personengruppen der gleichen Be- 
 schaftigungsart wurden in Amerika besonders wahrend des letzten Krieges einge- 
" fiihrt. Bei uns erschienen sie zuerst oft nur unter merkantilen Aspekten und daher 
“nicht immer iiberzeugend. Inzwischen sind auch hier die niitzlichen Seiten fiir 
~ unsere Verhdltnisse erkannt worden. Eine grobe Einteilung derartiger Messverfah- 
_ ren kann etwa folgende Dreiteilung aufweisen : 


1) Polyoptic AG., Basel. 
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A B Cc 
Messungen nach Schema mit Messungen durch optometrisch Messungen ~ 
angelerntem Experten, ohne geschulten Experten unter Be- durch 
besondere optometrische tiicksichtigung optischer Kor- Arzte 
Kenntnisse rekturmoglichkeiten 


Bei allen drei Gruppen wird vorausgesetzt, dass der Experte mit den Arbeitsbedin- 
gungen des Priiflings weitgehend vertraut ist. 

Unsere bisherige Erfahrung lehrt, dass sich Messungen nach «A» hauptsachlich 
fiir selten vorkommende Massenuntersuchungen (Militar) mit groben Ausscheidun- 
gen eignen und dass Priifungen durch Arzte kaum durchfiihrbar sind. Somit 
werden Priifungen mit optometrisch geschulten Experten «B», in kleineren Gruppen 
und mit etwas grésserem Zeitaufwand, am aussichtsreichsten. Es ist auch nicht zu 
iibersehen, dass die iibliche arztliche Betreuung unserer Bevélkerung schon eine 
Auswahl herbeifiihrt, welche grobe Ausscheidungsverfahren ohnehin hinter sich 
lasst. 

Wenn andererseits die Méglichkeit besteht, die Ergebnisse der Ortho-Rater- 
Priifungen «B» mit einem Arzt durchzugehen, so lasst sich die gesamte Basis ent- 
sprechend erweitern. 

Der Ortho-Rater ist iiblicherweise fiir folgende Messverfahren ausgeriistet: 
Messungen der Abweichungen der Augenachsen, unter Beriicksichtigung der Fu- 
sionsanregung. Messung der Sehscharfen mit einer Doppeltestmarke (Figur), sowohl 


binokular wie auch monokular bei binokularer Betrachtung. Ausserdem sind 
Stereopriifungen und Farbpriifungen méglich. Besonders diese letzteren verlangen 
vom Experten sichere Einfiihlung in die Sehbedingungen des Priiflings und eine 
gute Beherrschung der Materie, um gegebenenfalls feinere Methoden zu beriick- 
sichtigen. 

Weiter sind fiir geschulte Experten Refraktionsmessungen und einfache Prii- 
fungen des Gesichtsfeldes durchfiihrbar. In vielen Fallen werden vor allem die 
Refraktionspriifungen iiberhaupt nicht zu umgehen sein. Alle Messungen werden 
gegebenenfalls mit schon vorhandenen Brillen durchgefiihrt, wobei Riickschliisse 
auf Glaserzentrierung ins Gewicht fallen (Bifokalglaser). Aus dem gleichen Grunde 


werden auch die Testmarken fiir Ferne und Nahe mit unterschiedlichen Neigungs- 
winkeln der Blickrichtungen prasentiert. 5 
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Der bahnarztliche Dienst der SBB unter Leitung von Dr. SERATI hat eines der 
Ortho-Rater-Instrumente zusatzlich noch fiir spezifische Aufgaben dieses Gebietes 
vervollstandigen lassen. Der Apparat gestattet die Einfithrung farbiger Stereo- 
aufnahmen von Bahn- und Signalanlagen mit eingesetzten Lichtern, sowie Mes- 
sungen farbiger Signalbilder in variablen Gréssen und Kombinationen. 

Es zeigt sich, dass die Einfiihrung solcher verfeinerter Priifbedingungen den 
Fragenkomplex teilweise in neuem Lichte erscheinen lasst und dass unsere bisheri- 
gen Erkenntnisse iiber die erforderliche Erfahrung des Experten und die daraus 
abzuleitenden optometrischen Voraussetzungen bestatigt werden. 

_ Eine Ausdehnung solcher Messungen auf die Interessengruppen des Strassen- 
verkehrs scheint ebenfalls aussichtsreich, doch gelten auch dort sinngemass die 
angedeuteten Moglichkeiten und Grenzen. 


Ein GerAt zur serienmAssigen Priifung des Farbstichs von Objektiven’). 
Von W. Lotmar und R. Want, Aarau?) 
_ Moderne Photoobjektive zeigen haufig einen leichten spektralen Gang der 
Transmission, was einerseits von Absorption der verwendeten Glaser, andererseits 
von der Wellenlangenabhangigkeit der Wirkung aufgebrachter Reflexverminde- 
rungsschichten herriihrt*). Bei Farbaufnahmen kann dadurch eine Veraénderung 
der Tonwerte hervorgerufen werden, was vor allem dann auffallt, wenn die zu 
einer Kamera gehérigen Wechselobjektive in dieser Hinsicht differieren. 
_ Die Transmissionskurven von Objektiven, bei denen einer oder beide der er- 
wahnten Einfliisse wirksam sind, zeigen in der Regel den in Figur 1 wiedergege- 
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Figur 1 
Transmission eines Photoobjektivs, schematisch. 


pbenen Typus, das heisst die Anderung der Transmission ist im ganzen Bereich 
relativ gering. Zur Charakterisierung der photographischen Wirkung geniigt es 
daher, wie die Erfahrung gezeigt hat, das Verhaltnis zweier Werte im roten und 
blauen Spektralgebiet zu kennen, wobei die Auswahl der Wellenlangen nicht 


| 2) Ein ahnliches Gerat zur Produktionsiiberwachung von Interferenzschichten in Reflexion 
wurde von A. C. TRauB, JOSA 46, 999 (1956) beschrieben. 
8) Kern & Co, AG., Aarau. 

4) Siehe z. B. A. E. Murray, JOSA 46, 790 (1956). 
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kritisch ist. Die Transmissionswerte der durch die beiden Schottfilter RG 1 und. 
BG 12 (Figur 2) begrenzten Spektralbereiche erwiesen sich als hierfiir geeignet. 
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Figur 2 
Transmission der Schottfilter BG 12 und RG 1. 


Es wurde ein Gerat gebaut, welches die Differenz der Transmissionswerte eines 
Priiflings fiir Rot und Blau direkt und sofort anzeigt. Um bei Priiflingen verschie- 
dener Art vergleichbare Werte zu erhalten, wurde die Anlage mit einer automati- 
schen Pegelregulierung fiir das Rotsignal ausgestattet. 

Der Aufbau ist folgender (Figur 3): 


Automatische 
Pegelstabilisierung 


(Verstarkungsgrad- 
regelung) 


Amplituden - 
Vergleichs - 
Schaltung 


Spitzenwert- 
rot Gleichrichter 


lampen- = 
Sfabilisator 
~ / = 
elektronischer Verstarker Anzeige 
rotierende = Ly Umschalter 
Farbfilferscheibe = 


Synchronmotor . 
(3000 U/min) 


Figur 3 
Schema des Gerates. 


_ Der von einer stabilisierten 9-W-Gliithlampe ausgehende Lichtstrom durchsetzt 
einen Kondensor, eine rotierende F ilterscheibe mit einem Rot- und einem Blau- 


: 
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filter, hierauf den Priifling, und fallt auf ein Selen-Photoelement. Dessen Strom- 
stésse werden verstarkt und durch einen von derselben rotierenden Scheibe licht- 
elektrisch gesteuerten Umschalter in ein Rot- und ein Blausignal getrennt, deren 
Amplitudendifferenz dann tiber eine Vergleichsschaltung angezeigt wird. 

_ Als Kondensor wurde ein farbkorrigiertes Projektionsobjektiv verwendet. Vor 
das Photoelement wurde eine beidseitig mattierte Linse geschaltet, um eine még- 
lichst gleichmassige Lichtverteilung tiber die ganze Flache zu erreichen. 

Die Anzeigezeit des Geradtes (Intervall vom Einbringen des Priiflings bis zum 
Einspielen der Galvanometernadel) betragt 1,5 sec, die Empfindlichkeit 2,5 Skalen- 
teile fiir eine Anderung des Blausignals um 1%. Die Anzeige des Farbstichs kann 
geregelt werden und wird bei Abwesenheit eines Priiflings auf die Mitte der Skala 
(0) eingestellt. Als Standard zur Eichung des Ausschlages dient eine aus gelbstichi- 
gem Glas hergestellte Planplatte. 


Summary 


___ A photometer evaluating the transmission ratio of optical systems for two spectral 
‘ranges (red and blue) is described. Light falling alternately through two suitably 
‘chosen rotating spectral filters impinges on a selenium photoelement. Two inter- 
laced pulse trains are thus generated, which, after amplification, are electronically 
‘separated and compared. Sensitivity is 2-5 scale divisions for a ratio variation of 1%. 


_ Gesichtspunkte bei der Stoffauswahl fur dunne Schichten in der Optik. 
Von E. Ritter, Balzers®), 

_ Mit Hilfe absorptionsfreier diinner Schichten mit bestimmtem Brechungsindex 
und geeigneter Dicke lassen sich Interferenzeffekte erzielen, die in der Optik in 
“gunehmendem Masse verwendet werden (zum Beispiel Reflexionsverminderung, 
‘Lichtteilung, Kaltlichtspiegel, Schmalbandinterferenzfilter). Den theoretischen Ent- 
wiirfen stehen jedoch in der praktischen Ausfithrung zahlreiche Schwierigkeiten 
im Wege, besonders infolge der beschrankten Zahl brauchbarer Stoffe. Die Be- 
‘schrankung der Stoffzahl ergibt sich aus der Beriicksichtigung der folgenden 
_Gesichtspunkte: 

1. Herstellbarkeit des Stoffes in Form diinner Schichten, vor allem durch Ver- 
-dampfung im Hochvakuum. Etwa 130 Elemente und anorganische Verbindungen 
_(ausgenommen Legierungen) wurden bisher als diinne Schichten mittels Vakuum- 
-aufdampfung hergestellt. 

: Relativ leicht aufdampfbare Stoffgruppen sind: Elemente, Halogenide, Oxyde 

und Sulfide. 

- 2. Bestandigkeit gegen mechanische Beanspruchung und Haftfestigkeit. Be- 
sonders kritisch bei dicken Schichtpaketen fiir UR-Systeme. Durch Heizung des 

_Tragers wahrend der Bedampfung giinstige Beeinflussung. 

. 3. Bestandigkeit gegen chemische Angriffe, vor allem gegen Feuchtigkeit. 

“Diese Forderung beschrankt besonders die im UV verwendbaren Stoffe. 

4, Gewiinschte optische Eigenschaften. 

__a) Absorptionsfreiheit, das heisst bei diinnen Schichten fiir den sichtbaren Be- 

Breich K < 10? cm". . 

_ _b) Méglichst konstanter Brechungsindex iiber méglichst weiten Bereich. Untere 

Grenze stets 1,30-1,34. Obere Grenze im UV 2,05-2,1, im Sichtbaren 2,4, im UR 

3,5—-5,4. 


5) Balzers Aktiengesellschaft fiir Hochvakuumtechnik und diinne Schichten, Balzers, Fiirsten- 
~ tum Lichtenstein. 
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Die nachfolgende Tabelle fasst die Stoffe zusammen, die sich bisher bewahrt 


haben. Die angegebenen Zahlen sind lediglich Richtwerte: 


Stoff n(A) Absorptionsfreier 
Bereich | 

Na,AlF, 1,33 (550 my) 250 my—14 pu 

MgF, 1,38 (550 mp) 250 my-10 

Si,O,-SiO, 1,55 (550 mp) 300 my— 8 u 

SiO 1,6-1,9 (3 1) 800 mu— 8 4 

ThF, 1,52 (550 mp) 250 mu-12 pu 

LaF, 1,75 (350 my) 1,6 (550 mp) ab 250 mu 

NdF, 1,75 (350 mp) 1,6 (550 my) ab 250 mp 

CeF, 1,75 (350 mp) 1,6 (550 mp) ab 250 mu 

La,O; 2,0 (400 my) 1,85 (550 mp) ab 300 mu 

Nd,O, 2,1 (400 mp) 2,0 (550 mp) ab 300 mu 

PbF, 2,05 (250 mu) 1,75 (550 mp) 250 my-17 

ceo, 2,2 (550 mp) 400 my-16 

ZnS 2,3, (550 my) 350 my—14,5 wu 

AMO} 2,4 (550 mp) 390 myu-12 wu 

Si 3,5 (1,4 p) ab 800 mu 

Ge ale (24) abr isi 

Te 5,4 (3,5 p) ab 3,5) 


175 Jahre Verlag Friedr. Vieweg & Sohn, Braunschweig 


Friedrich Vieweg griindete am 1. April 1786 den Verlag, der noch heute seinen 
Namen tragt. Zu seinen Autoren zahlen unter anderen HERMANN Vv. HELMHOLTZ, 
Madame Curiz, ARNOLD SOMMERFELD, Lupwic PRANDTL, ALBERT EINSTEIN, 
Max v. Lave und ERwIN scHRODINGER. Aus Anlass dieses Jubilaums stiftet der 


Verlag Preise von insgesamt DM 15 000.— fiir wissenschaftliche Arbeiten auf den | 
Gebieten der Mathematik, Physik und Chemie, die sich zur Ver6ffentlichung in — 


Buchform eignen. Dem Kuratorium dieser Stiftung gehéren an: die Professoren 
J. BarTELs, W.GeERLacH, W. Haack, R. Huiscen, J. Mattraucn, W. QUADE, 
F, Sauter, F. See, H. SIEDENTOPF und E. WICcKE. Die Ausschreibungsbedingun- 
gen stellt der Verlag Interessenten auf Wunsch gern zur Verfiigung. 


IFIPS-Kongress, Miinchen 1962 


Die International Federation of Information Processing Societies wird vom 
27. August bis zum 1. September 1962 in Miinchen einen Kongress abhalten. Vor- 
tragsanmeldungen, begleitet von einer Zusammenfassung in englischer Sprache im 
Umfange von 500-1000 Worten, sind bis 15. September 1961 einzureichen, entweder 
an den Obmann des Programmausschusses, Herrn Prof, Dr. F. L. BaveEr, Institut 
fiir Angewandte Mathematik der Johannes-Gutenberg-Universitat, Mainz, Jakob- 
Welder-Weg 7; oder an den Stellvertreter, Herrn Dr. H. BILiine, Max-Planck- 
Institut fiir Physik und Astrophysik, Institut fiir Astrophysik, Abt. Numerische 
Rechenmaschinen, Miinchen 23, Aumeisterstrasse. Die Vortragsanmeldungen wer- 
den auf Grund der Zusammenfassungen vom internationalen Programm-Komitee 
der IFIPS gepriift und alsdann werden die Autoren der angenommenen Zusammen- 
fassungen aufgefordert, vollstandige Manuskripte vorzulegen. 


i 
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Asymptotic Behavior and Stability Problems in Ordinary Differential 
Equations. Von LAMBERTO CrEsaArRI. Ergebnisse der Mathematik und ihrer Grenz- 
gebiete. Neue Folge, Heft 16 (Springer-Verlag, Berlin-Géttingen-Heidelberg, 1959). 
271 S., 37 Fig.; DM 68.-. 

Von der Theorie der gewohnlichen Differentialgleichungen sind fiir die Anwen- 
dungen Aussagen iiber das asymptotische Verhalten und die Stabilitat der L6- 
sungen wesentlich wichtiger als die den reinen Mathematiker interessierenden 
Fragen der Existenz und Eindeutigkeit. Die modernen Entwicklungen der Technik, 
‘speziell der Regelungstechnik und der Automation, haben in den letzten Jahr- 
zehnten zu einer Fiille von diesbeziiglichen Untersuchungen Anlass gegeben, deren 
iibersichtliche Zusammenfassung und Darstellung keine leichte Aufgabe ist. Dem 
Bericht von L. Crsari ist die Lésung dieser Aufgabe aber in vortrefflicher Weise 
gelungen. Im ersten Kapitel werden die verschiedenen Stabilitatsbegriffe diskutiert 
und die linearen Systeme mit konstanten Koeffizienten betrachtet. Das folgende 
‘Kapitel behandelt die allgemeinen linearen Systeme und bringt die Theorie von 
Lyapunov fiir Systeme mit variablen Koeffizienten, die Theorie von FLogueEt fir 
solche mit periodischen Koeffizienten und ausfiihrlich die Gleichungen 2. Ordnung. 
Das dritte Kapitel, das nahezu die Halfte des Buches ausmacht, ist den nicht- 
linearen Systemen gewidmet und entwickelt die erste und zweite Methode von 
Lyapunov, die analytischen Methoden (Porncaré, KryLoyv-BocoLjuBov, Sto- 
rungsmethoden) und endlich die topologischen Methoden (PoINcaRE, BENDIXSON). 
Das abschliessende Kapitel bringt dann noch die asymptotischen Entwicklungen. 
Der Bericht hebt die allgemeinen Linien hervor und bringt in jedem Abschnitt 
neben den allgemeinen Theoremen meist nur einige wenige typische Beweise, dazu 
illustrative Beispiele. Das bis 1957 reichende Literaturverzeichnis beansprucht nicht 
weniger als 70 Seiten. Das Buch vermittelt jedem, der sich mit diesen Fragen zu 
befassen hat, einen ausgezeichneten Uberblick und die entsprechenden Hinweise 
und wird ihm ein unentbehrlicher Helfer sein. E. RotH-DESMEULES 


Partielle Differentialgleichungen. Von G. Hettwie (B. G. Teubner, Ver- 
lagsgesellschaft, Stuttgart 1960). 346 S., 35 Fig.; DM 29.80. 

Im Rahmen der von G. Kérue betreuten «mathematischen Leitfaden» des 
Teubner-Verlags gibt der Verfasser eine Einfiihrung in die partiellen Differential- 
gleichungen, die vor allem den theoretisch-mathematischen Standpunkt betont, 
also sich mit Existenz- und Eindeutigkeitsfragen befasst. Dabei hat er sich bemiuht, 
auch die neuesten Forschungen zu Wort kommen zu lassen, was ihm in didaktisch 
hhervorragender Weise gegliickt ist. Nach einem Rundblick iiber die in der mathe- 
‘matischen Physik auftretenden partiellen Gleichungen wird ausfiihrlich die Ein- 
teilung in elliptische und hyperbolische Gleichungen besprochen, die Charakteri- 
‘stiken werden eingefiihrt und Normalformen hergeleitet. Zur Abklarung der Ein- 
deutigkeitsfragen wird in origineller Weise hauptsdchlich die Energieintegral- 
‘methode herangezogen. Sehr modern ist das Kapitel tiber Existenzfragen geschrie- 
‘ben. Die Anfangsprobleme beim hyperbolischen Typus werden mittels der Charak- 
teristikentheorie erledigt, wahrend fiir die gemischten Probleme die Laplace- 
Transformation herangezogen wird. Bei den elliptischen Randwertproblemen folgt 
‘der Verfasser den neueren Methoden von FRIEDRICHS, Lax, NIRENBERG, die Ver- 
feinerungen der direkten Methoden der Variationsrechnung sind. Ausserdem sind 
einige speziellere Verfahren geschildert (Methode der spharischen Mittelwerte, un- 
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publizierte Untersuchungen von WIENHOLTz tiber die Annahme der Randwerte, 
Weylsches Lemma). Ein letztes Kapitel bringt eine Einfiihrung in die eh | 
analytischen Methoden und in die Schaudersche Beweistechnik. ee | 
Es ist ein grosses Verdienst des Verfassers, die verstreute neuere Literatur 
iiber partielle Differentialgleichungen iiberarbeitet und in einem angenehm lesbaren 
Buch zusammengestellt zu haben. Das Werk wird auch dem angewandten Mathe- 
matiker willkommen sein, der sich um eine theoretische Klarung seiner Probleme. 


bemiiht. _ E. STreFEL 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLuccr, Band 37/1: Atome III/Molekiile I (Springer-Verlag, Berlin 1959)§ 
439 S., 215 Fig.; DM 96.-. 

Der vorliegende Teilband Atome III — Molekeln I der Gruppe VII des Hand-. 
buchs (Atom- und Molekiilphysik) enthalt nur drei Kapitel. Von diesen beansprucht 
das erste iiber Atom- und Molekularstrahl-Spektroskopie (P. Kuscu und V. Wa} 
HucuEs etwa 170 Seiten) zurzeit hauptsachlich das Interesse des Physikers, wah- 
rend die beiden andern iiber die Schwingungs- und Rotationsenergie von Molekiilen | 
(von H. H. Nrersen, etwa 110 Seiten) bzw. iiber Elektronenstossprozesse mit 
Molekeln (von J. D. Craccs und H. S. W. Massey) eher der heutigen physikali- 
schen Chemie zugeordnet werden kénnen. 

Der Beitrag tiber Molekularstrahl-Spektroskopie ist ahnlich, wenngleich be- 
deutend konzentrierter, dem bekannten Buch von N. F. Ramsey und in drei 
Abschnitte iiber experimentelle Methodik, Theorie und Ergebnisse gefasst. Dariiber 
hinaus ist die Literatur bis und mit 1958 beriicksichtigt und die Theorie ist starker | 
betont. Man darf diesen Beitrag als hervorragende Ubersicht bezeichnen, fiir 
einen Handbuchartikel scheint er recht knapp konzipiert. 

H. H. NreLsen’s Beitrag iiber Rotations- und Schwingungsenergie von Mole- 
keln ist in dieser Form unseres Wissens erstmalig. Die zahlreichen Arbeiten iiber 
die Quantenmechanik der Rotations-Schwingungszustande von Molekeln — eine 
grosse Anzahl stammt vom Autor selbst — sind hier unter einheitlichem Gesichts- 
punkt zusammengefasst. Ausgehend von der klassischen Hamiltonfunktion von 
rotierenden Massenpunktsystemen mit kleinen Schwingungen, aber ohne interne 
Rotationen, wird deren Quantenmechanik nach steigenden Naherungen speziell 
der Wechselwirkung zwischen Rotation und Schwingungsenergie und der Anhar- 
monizitat behandelt. Dabei werden Rotationsprobleme des starren Kreisels und 
das Schwingungsproblem in harmonischer Naherung formal behandelt, ohne Wie- 
dergabe der bekannten numerischen Tabellen der Energieeigenwerte des asymme- 
trischen Kreisels. Fiir die Fille der symmetrischen Kreisel mit Schwingungsdreh- 
impuls sowie fiir die héheren N daherungen des Schwingungsrotationsproblems 
werden (in der Nielsenschen Formulierung) sehr niitzliche Tabellen fiir die Ele- 
mente der Hamiltonmatrix gegeben. Ferner wird die Anwendung der Theorie auf 
die Analyse von Infrarotspektren hoher Auflésung eingehend demonstriert, speziell 
die praktisch oft beobachteten Anomalien, die durch L-type doubling, anharmo- 
nische und rotatorische Resonanz (bis zur 2, Naherung) verursacht werden. 

Auch der dritte Beitrag iiber Elektronenstossprozesse mit Molekeln erscheint 
in so kompakter Form wohl zum ersten Mal in der Literatur. Die (knappe) Darstel- 
lung der Theorie der Stossprozesse zwischen schnellen und langsamen Elektronen 
und Molekeln und die dadurch verursachten Ubergange der Molekelzustande 
findet unmittelbare Anwendung auf eine grosse Zahl experimenteller Daten. 
Dadurch besitzt der Artikel fiir den Massenspektroskopiker nicht nur kompilatori- 
schen Charakter, sondern er kann geradezu als Anleitung zur Analyse der Massen- 
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spektren von Molekeln betrachtet werden. Gemessen am Stand der Kenntnisse auf 
diesem Gebiet kommt die Konzeption dieses Beitrags dem Ideal eines Handbuch- 
artikels gewiss sehr nahe. Hs. H. GUNTHARD 


“ 


_ Confluent Hypergeometric Functions. Von L. J. SLaTER (Cambridge Uni- 
versity Press, London 1960). 247 S., 18 Abb.; 65s. 

In der Physik wird man gelegentlich auf die konfluente hypergeometrische 
Funktion gefiihrt, wenn man versucht, eine partielle Differentialgleichung durch 
Separation in Kugelkoordinaten zu lésen. Ein bekanntes Beispiel ist die quanten- 
mechanische Behandlung des Wasserstoffatoms. Die konfluente Funktion ist einer- 
seits Grenzfall der allgemeinen GauS8schen hypergeometrischen Funktion und ent- 
halt andererseits Besselfunktionen, Whittakersche Funktionen und Laguerresche 
Polynome als Spezialfalle. 

Die vorliegende Monographie leitet zunachst die Funktionaleigenschaften der 
konfluenten Funktion her, wobei von der grundlegenden Kummerschen Differen- 
tialgleichung ausgegangen wird. Die beiden Fundamentallésungen dieser linearen — 
und im Nullpunkt singuléren — Differentialgleichung 2. Ordnung werden sorgfaltig 
diskutiert und die Beziehungen zu den Whittakerfunktionen hergestellt. Ein zwei- 
tes Kapitel ist den Differentialeigenschaften und Rekursionsformeln gewidmet, die 
konfluente Funktionen F(a, b; x) verkniipfen, die zu verschiedenen Werten der 
Parameter a, b gehoren. Ein drittes Kapitel bringt die wichtigen Darstellungen der 
konfluenten Funktion durch Umlaufintegrale in der komplexen Ebene, einschliess- 
lich der Laplace-Transformierten von F(a, b; x). Diese Integraldarstellungen bilden 
das Instrument fiir die asymptotischen Entwicklungen der konfluenten Funktion, 
die im 4. Kapitel des Werkes in grosser Vollstandigkeit auseinandergesetzt sind. 

Besonders ansprechend ist der originelle Versuch des Verfassers, in einem 
Abschnitt «descriptive properties», dem Leser die konfluente Funktion auch an- 
schaulich naherzubringen. Man findet dort graphische Darstellungen der Nullstel- 
Jenkonfigurationen und Reliefs des Verlaufs der Funktion im Komplexen. 

Speziell sei noch auf die Tabellen hingewiesen, die dem Buch beigegeben sind. 
Numerische Tabellen von F(a, b; x) miissen naturgemass rudimentar bleiben, da 
die Funktion von 3 Variablen a, b, x abhangt. Man findet aber mehr Information 
als bisher in der Literatur, so zum Beispiel eine Tabelle von F(a, b; x) mit dem 
Schritt 0,1 in allen drei Variabeln. a lauft von (— 1) bis + 1, 6 von 0,1 bis 1 und 
x von 0,1 bis 10. Eine weitere Tabelle der speziellen Werte F(a, b; 1) weitet altere 
Tabellen von MILLER aus. 

Das Werk von SLATER erganzt den Ergebnisbericht von H. BucnHo1z (Die 
honfluente hypergeometrische Funktion, Springer 1953), der viel benutzt wird, und 
4xann niitzliche Dienste neben der bekannten Formelsammlung von ERDELYI, 
Macnus, OBERHETTINGER, TRiIcom1 (Higher transcendental functions, Mc Graw- 
‘Hill, 1953) leisten. E. STIEFEL 


Die Mechanisierung des Weltbildes. Von E. J. D1)KSTERHUIS. Ubersetzt 
aus dem Hollandischen von HELrca HaBIcHT (Springer-Verlag Berlin, Gottingen, 
Heidelberg 1956). 594 S., 47 Abb.; DM 36.—. 

Auch wenn das hier anzuzeigende Werk von DIJKSTERHUIS nicht direkt zu den 
Gebieten der ZAMP gehort, so diirfte es doch fiir manchen Leser von erheblichem 
Interesse sein. Gerade in der heutigen Zeit, wo man — wenn es so ausgedriickt 
werden darf — ohne grosse Hemmungen neue Theorien aufstellt und diese allein 
‘nach ihrer Brauchbarkeit und nach ihrem Erfolg, die Phanomene darzustellen und 
-vorauszusagen, beurteilt werden, scheint ein Riickblick auf die ausserordentlich 
Jange Entwicklungszeit der modernen Naturwissenschaften von grosser Bedeutung 
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zu sein, D1jKSTERHUIS’ Darstellung — auf deren ungeheure Fiille hier nicht einm: : 
andeutungsweise eingegangen werden kann — beginnt mit den griechischen Denkern, 
den ersten Begriindern wissenschaftlicher Theorien und endet mit Newtons Philo- 
sophiae Naturalis Principia Mathematica. In meisterhaft klarer und knapp@ 
Weise werden die Hauptstrémungen der griechischen Philosophie dargestellt, 
sowie die verschiedenen Fachwissenschaften wie Mathematik, Astronomie, Physik, 
Technik, Chemie, Astrologie, die die Grundlage bilden fiir die spatern Entwicklun- 
gen. Nach einer langen, mehr oder weniger unfruchtbaren Ubergangsperiode 
beginnt im 12. Jahrhundert ein Wiedererwachen des naturwissenschaftlichen Inter- 
esses in Europa. Durch THomas von AguiN erfolgte dann die Synthese der christ- 
lichen Lehre mit der Philosophie von ARISTOTELES, wodnrch die Theologie praktisch 
die Oberaufsicht tiber alle Wissensgebiete erhielt, was dann spditer einerseits zu 
zahlreichen und wesentlichen Konflikten zwischen Kirche un? Wissenschaft fiihren 
musste und andrerseits dem Denken des Stagiriten eine beherrschende Rolle zuwies, 
von der man sich nur allmahlich und mit grésster Miihe wieder lésen konnte. Hu- 
manismus und Renaissance bereiten dann die klassische Naturwissenschaft vor, die 
mit KoPERNIKUs etwa beginnt und iiber KEPLER, GALILEI, HuyGENsS, DESCARTES 
schliesslich zu Isaac Newton hinfiihrt. In dieser Zeit erhalten Hydrostatik, Optik: 
und Mechanik allmahlich eine mathematische Gestalt und die empirische Methode: 
gewinnt die ihr zukommende Bedeutung, womit dann der Schritt zur klassischen. 
und damit auch zur modernen Naturwissenschaft im Prinzip vollzogen wurde. 

Das Buch von DijKsTERHUIS zeigt in eindriicklicher und fesselnder Weise die. 
vielen Um- und Abwege des menschlichen Denkens, bis endlich die Befreiung von 
vorgefassten philosophischen und religidsen Ansichten méglich wurde und der 
geeignete Weg zur Erforschung der Natur und ihrer Gesetzmassigkeiten gefunden 
wurde. Wesentlich ist — was eigentlich selbstverstandlich sein sollte, aber doch 
besonders betont werden muss — dass man wieder gelernt hat, zu den Quellen 
zurtickzukehren, um das historische Studium hier zu beginnen; denn wie manches. 
wurde im Laufe der Zeit falsch oder von einem ungeeigneten Gesichtspunkt aus inter- 
pretiert und wie schwierig ist es auch heute noch 6fters, die tatsdchliche Meinung 
alterer Autoren zu ergriinden (wie manche Kontroversen zeigen). 

DijKsTERHUIs, der sich selbstverstandlich auf zahlreiche andere Arbeiten aus 
der neueren Zeit stiitzen muss, raumt mit vielen irrigen Ansichten auf, die noch in 
zahllosen Lehrbiichern der Physik oder der Naturwissenschaften zu finden sind 
von Verfassern, die sich verpflichtet fiihlen, historische Bemerkungen einzuflechten 
und sich dabei auf Darstellungen aus n-ter Hand (mit n > 4) stiitzen. Im ganzen 
gesehen kommt das Mittelalter dabei doch bedeutend besser weg, als die landlaufige 
Meinung es wahrhaben will; und vieles auf dem Gebiete der Botanik, Zoologie und 
Medizin ist aus den vorhandenen Handschriften — von denen moderne kritische 
Ausgaben weitgehend fehlen — noch gar nicht erschlossen worden, und sie diirften 
das Bild des «dunklen» Mittelalters noch wesentlich erhellen. — Man muss DijKs- 
TERHUIS dankbar sein, dass er die gewaltige Miihe auf sich genommen hatte, die 
Ergebnisse der modernen historischen Forschung auf dem Gebiete der exakten 
Naturwissenschaften in einer grossartigen Uberschau in ihren wesentlichen Ziigen 
zusammenzufassen — und dies ohne den Leser durch seitenlange Zitate und Anmer- 
kungen zu langweilen, und trotzdem seiner Pflicht der Quellenangabe in ausrei- 
chendem Masse nachkommend. Das Studium dieses ausserordentlich griindlichen 
Werkes lohnt sich fiir jeden nicht ausschliesslich fachwissenschaftlich orientierten 
Naturwissenschaftler, und im Hinblick auf die heute geradezu tragische Bedeutung 
und Konsequenz mancher wissenschaftlicher F orschung fiir die Menschheit ist auch 
eine solche historische Besinnung geradezu eine Notwendigkeit. 
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INHALT: Kapitel I — Unendliche Mengen. Allgemeine Mengen — Punktmengen — Weitere 
Diskussion der infinitaren Eigenschaften von Mengen. Kapitel II — Funktionen auf M engen, 
Konvergenz von Funktionen auf Punktmengen — Stetigkeit von Funktionen auf Punkt- 
mengen. Kapitel III — Unendliche Folgen und Rethen. Unendliche Folgen — Unendliche 
Reihen — Funktionenfolgen und Funktionenreihen — Potenzreihen. Kapitel IV — Ergan- 
zungen zur Differentialrechnung. Weitere Anwendungen der Differentialrechnung auf Diskus- 
sion der Funktion einer Variablen — Differentiation bei Funktionen mehrerer Variablen — 
Partielle Ableitungen héherer Ordnung — Jacobische Matrizen und Determinanten — Das 
totale Differential. Kapitel V — Anwendungen der Differentialrechnung auf die Analysis. 
Existenz der Lésungen von Gleichungen und Gleichungssystemen. Implizite Funktionen — 
Umkehrung von Abbildungen und Funktionensystemen — Extrema bei Funktionen mehrerer 
Variablen. Kapitel VI — Numerische Rechenmethoden. Anwendungen der Differentialrechnung 
auf Interpolation — Numerische Differentiation und Integration — Angenaherte Auflésung 
von Gleichungen — Bernoullische Zahlen und Polynome — Die Euler-Maclaurinsche Formel. 
Kapitel VII — Allgemeines tiber Kurven. Ebene Kurven. Grossen erster Ordnung in der 
Kurventheorie — Gréssen zweiter Ordnung in der Kurventheorie — Evolute, Evolvente und 
Parallelkurven — Enveloppen von Kurvenscharen. Kapitel VIII — Raumkurven und Flachen. 
Vektoren — Die Schmiegungsebene einer Raumkurve — Kriimmung und Torsion — Tangen- 
tialebene und Flachennormale — Weitere Ausfiihrungen tber Flachen (Enveloppen; Geo- 
metrie auf der Flache). Register. 

In dieser Auflage ist der zweite Band stark entlastet worden, da die Aufgaben ganz 
herausgenommen wurden und einige Teile der Reihenlehre und der Kurventheorie bereits 
in der Neuauflage des ersten Bandes behandelt worden sind. So ergab sich die Méglichkeit, 
bei der Neubearbeitung verschiedene Dinge einzubeziehen, die die Darstellung der heu- 
tigen Entwicklung der Mathematik und ihrer Anwendungen besser anzupassen gestatten. 
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